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Abstract 

We provide a comprehensive, up-to-date analysis of possible New Physics contributions to the 
mass difference AMd in D°-D mixing. We consider the most general low energy effective Hamil- 
tonian and include leading order QCD running of effective operators. We then explore an extensive 
list of possible New Physics models that can generate these operators, which we organize as in- 
cluding Extra Fermions, Extra Gauge Bosons, Extra Scalars, Extra Space Dimensions and Extra 
Symmetries. For each model we place restrictions on the allowed parameter space using the recent 
evidence for observation of D meson mixing. In many scenarios, we find strong constraints that 
surpass those from other search techniques and provide an important test of flavor changing neutral 
currents in the up-quark sector. We also review the recent BaBar and Belle findings, and describe 
the current status of the Standard Model predictions of D^-D^ mixing. 



1 



I. INTRODUCTION 



Meson-antimeson mixing has traditionally been of importance because it is sensitive to 
heavy degrees of freedom that propagate in the underlying mixing amplitudes. Estimates 
of the charm quark and top quark mass scales were inferred from the observation of mixing 
in the and Bd systems, respectively, before these particles were discovered directly. 

This success has motivated attempts to indirectly detect New Physics (NP) signals by 
comparing the observed meson mixing with predictions of the Standard Model (SM). Mixing 
in the Kaon sector has historically placed stringent constraints on the parameter space 
of theories beyond the SM and provides an essential hurdle that must be passed in the 
construction of models with NP. However, anticipated breakthroughs from the B-factories 
and the Tevatron collider have not been borne out - the large mixing signal in the and Bg 
systems is successfully described in terms of the SM alone (although the parameter spaces 
of various NP models have become increasingly constrained). Short of awaiting LHCB and 
the construction of a super-B facility, there is one remaining example for possibly observing 
indirect signs of NP in meson mixing, the flavor oscillations. In this case, the SM 
mixing rate is sufficiently small that the NP component might be able to compete 
There has been a flurry of recent experimental activity regarding the detection of D^-D^ 
mixmg 

mm 

, which marks the first time Flavor Changing Neutral Currents (FCNC) 
have been observed in the charged +2/3 quark sector. With the potential window to discern 
large NP effects in the charm sector p, Q] and the anticipated improved accuracy for future 
mixing measurements, the motivation for a comprehensive up-to-date theoretical analysis of 
New Physics contributions to D meson mixing is compelling. 



A. Observation of Charm Mixing 

The heightened interest in D^-D^ mixing started with the almost simultaneous observa- 
tions by the BaBar [5] and Belle collaborations of nonzero mixing signals at about the 
per cent level, ^ 

?/d = (0.97 ±0.44 ±0.31) ■ 10^2 (BaBar), (1) 
^ Our definitions of the mixing parameters x^, J/d, Vu and y'^^^ are standard and are given in Eqs. (|lip . (|12p . 
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yh 



fJ^P) = (1.31 ±0.32 ±0.25) ■ 10" 



(Belle) 



(2) 



This was soon followed by the announcement by the Be^ 
ments from the Dalitz plot analyses of —>■ i^'^Tr^vr 



le collaboration of mixing measure- 

3, 



xd = (0.80 ±0.29 ±0.17) ■ 10" 



1/D = (0.33 ±0.24 ±0.15) ■ 10" 



(3) 



A pre 



iminary fit to the current database^ by the Heavy Flavor Averaging Group (HFAG) 
gives ^ 



7/D = (6.6±2.1)-10- 



(4) 



Since this paper addresses the issue of the mass splitting induced by mixing, our primary 
concern is with the signal for xd, seen here to be a 2.4 sigma effect. This is below the 
generally accepted threshold for "evidence" and is more in the nature of a "hint" . However, 
we note that a 2.4 sigma effect will automatically have a non-zero lower bound at 95% 
confidence-level. For the sake of reference, we cite the one-sigma window for the HFAG 
value of xb, 

5.4-10"^ < xd < 11.7- 10"^ (one — sigma window) , (5) 
or equivalently for AMd itself. 



8.7 ■ 10"'^ GeV < AMd < 1-9 ■ 10"'* GeV (one - sigma window) 



(6) 



Let us briefiy describe our strategy for dealing with the above HFAG values in light 
of both SM and NP contributions. We shall argue in Sect. HI that the SM predictions, 
although indeed compatible with the observed range of values for the D mixing parameters, 
contain significant hadronic uncertainties. Moreover, we do not know the relative phase 
between the SM contribution and that from any NP model, so that xd will lie between the 
extreme limiting cases of constructive and destructive interference. In addition, since the 
observation of D mixing is new, the measurements will fiuctuate with future refinements in 
the analyses and as more data is collected. To best deal with these realities, we will present 



2 An updated fit[8| gives the values xu = SAtfl ■ 10"^ , yu = (6.9 ± 2.1) • 10~^. These are essentially 
unchanged from the HFAG preliminary results given above and used in our analysis; the difference will 
not affect our numerical results. 
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our results by displaying a given NP prediction as a pure NP signal {i.e. as if there were no 
SM component) and for comparison, display curves of constant xd for the five values 

xd = 15.0 ■ 10"^ , 11.7-10-^, 8.0-10-^, 5.0-10-^, 3.0 ■ 10"^ . (7) 

This (approximately HFAG 2cr) range reveals the sensitivity of xb to variations in the 
underlying NP parameter space. We will then show the present constraints placed on the 
NP model parameter space, by assuming that the NP contribution cannot exceed the la 
upper bound onxo- This procedure mirrors that which is traditionally employed in obtaining 
bounds on NP from mixing. 



B. New Physics Possibilities 

D^-D^ mixing at the observed level is much larger than the quark-level ('short-distance') 
SM prediction [Oj but is in quahtative accord with hadron- level ('long-distance') SM expec- 
tations. However, because the latter are beset with hadronic uncertainties, it cannot be 
rigorously concluded that only SM physics is being detected. In this paper, we will consider 
a broad menu of NP possibilities. As the operation of the LHC looms near, the number 
of potentially viable NP models has never been greater. Our organizational approach to 
analyzing these is to address NP models with: 

1. Extra Fermions (Sect. IV) 

A: Fourth Generation 

B: Heavy Vector-like Quarks 

(1) : g = -1/3 Singlet Quarks 

(2) : Q = +2/3 Singlet Quarks 
C: Little Higgs Models 

2. Extra Gauge Bosons (Sect. V) 

A: Generic Z' Models 

B: Family Symmetries 

C: Left-Right Symmetric Model 
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D: Alternate Left- Right Models from Eq Theories 
E: Vector Leptoquark Bosons 

3. Extra Scalars (Sect. VI) 

A: Flavor Conserving Two-Higgs-Doublet Models 
B: Flavor Changing Neutral Higgs Models 
C: Scalar Leptoquark Bosons 
D: Higgsless Models 

4. Extra Space Dimensions (Sect. VII) 

A: Universal Extra Dimensions 
B: Split Fermion Models 
C: Warped Geometries 

5. Extra Symmetries (Sect. VIII) 

A: Minimal Supersymmetric Standard Model 
B: Quark-Squark Alignment Models 
C: Supersymmetry with R-Parity Violation 
D: Split Supersymmetry 

In the above, we have chosen to consider only supersymmetry in Sect. VIII due to its 
extensive literature and to cover other extended symmetries elsewhere in the paper. 

Any NP degree of freedom will generally be associated with a generic heavy mass scale 
M, at which the NP interaction will be most naturally described. At the scale rric of the 
charm mass, this description will have been modified by the effects of QCD. These should 
not be neglected, so we perform our NP analyses at one-loop level for the strong interactions. 
The theoretical background for this is presented in Sect. II. 

Finally, in order to place the NP discussion within its proper context, it makes sense to 
first review SM charm mixing. This is done in Sect. III. The remainder of the paper then 
amounts to considering charm mixing with lots of 'extras'. The paper concludes in Sect. IX 
with a summary of our findings. 
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C. Basic Formalism 



Let us first review some formal aspects of charm mixing. The mixing arises from | AC| = 2 
interactions that generate off-diagonal terms in the mass matrix for and D° mesons. The 
expansion of the off-diagonal terms in the neutral D mass matrix to second order in the weak 
interaction is 

where i^J^*^'^^ and iJ]^*"'^^ are the effective |AC| = 2 and |AC| = 1 hamiltonians. 

The off-diagonal mass-matrix terms induce mass eigenstates Di and D2 that are super- 
positions of the flavor eigenstates and D'^, 

Di =pD^ ± , (9) 

2 

where |pp + |gp = 1. The key quantities in mixing are the mass and width differences, 

AMd = Ml - M2 and AFd = Ti - , (10) 
or equivalently their dimensionless equivalents, 

_amd , _ArD 

xd = — , and t/d = — ^ , (11) 

i D D 

where Fd is the average width of the two neutral D meson mass eigenstates. Two quantities, 
y^^^ and i/^, which are actually measured in most experimental determinations of AF^, are 
defined as 

y^^^ = (F+ - F_)/ (F+ + F_) = ?/D cos - xb sin - Aprod^ , 

y'l) = Vd cos 6Kn - xb sin Sktt , (12) 

where the transition rates r± pertain to decay into final states of definite CP, Ap^od = 
(^Njjo — Nj^^ / (^Njjo + A^^o^ is the so-called production asymmetry of and (giving the 
relative weight of and in the sample) and 6ktt is the strong pha se difference between 
the Cabibbo favored and double Cabibbo suppressed amphtudes [10|. The quantities Am 
and account for the presence of CP violation in D^-D^ mixing, with Am being related to 
the q,p parameters of Eq. ([9]) as Am = Iq/pI"^ — 1 and a CP-violating phase of M21 (if one 



neglects direct CP- violation) [ll|. In practice, ycp is measured by comparing decays of 



into a state of definite CP, such as K^K~ , to decays of into a final state which is not a 
CP-eigenstate (such as Kit) whereas y' is extracted from a time-dependent analysis of the 
D Kit transition fll]. 

The states Di allow for effects of CP violation. However, CP violation in mixing is 

n n n 

negligible in the Standard Model and there is no evidence for it experimentally [2|, Il2|, |13 |. 



Many New Physics scenarios contain new phases which can induce sizable CP violation in 
the D meson sector. Nonetheless, a thorough investigation of such effects is beyond the 
scope of the present paper. Therefore, we shall work in the limit of CP invariance (so that 
p = q) for the remainder of this paper. Throughout, our phase convention will be 

CVD° = + ^ . (13) 

Then Di become the CP eigenstates D± with CV D± = ± D±. 

2 

Keeping in mind the neglect of CP- violation and also the phase convention of Eq. (fT3ll . 
we relate the mixing quantities xd and ?/d to the mixing matrix as 



xr, = Re 

2MDrD 



2{W\H\''^\='\D') + (D°|z|d^a;T{7^L^^I=i(x)7^L^^I=i(0)}|D°) 



VB = Im J d'x T{nl^''\=\x) nl^^\=\0)}\D') , (14) 

where 7il^^l=i (x) is the weak hamiltonian density for |AC| = 1 transitions and T denotes 
the time-ordered product. There is no contribution to ?/d from the local |AC| = 2 term, as 
it has no absorptive part. New Physics contributions to ?/d have already been addressed in 
Ref. [14], so the primary thrust of this paper will be to focus on xd- 

The next step, in Sect. II, is to expand the time-ordered product of Eq. (fT^ in local 
operators of increasing dimension (higher dimension operators being suppressed by powers 
of AqcdMc)- 



II. GENERIC OPERATOR ANALYSIS OF D^-D^ MIXING 

Though the particles present in models with New Physics may not be produced in charm 
quark decays, their effects can nonetheless be seen in the form of effective operators generated 
by the exchanges of these new particles. Even without specifying the form of these new 
interactions, we know that their effect is to introduce several |AC| = 2 effective operators 
built out of the SM degrees of freedom. 
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A. Operator Product Expansion and Renormalization Group 



By integrating out new degrees of freedom associated with new interactions at a scale 
M, we are left with an effective hamiltonian written in the form of a series of operators 
of increasing dimension. Operator power counting then tells us the most important con- 
tributions are given by the operators of the lowest possible dimension, c? = 6 in this case. 
This means that they must contain only quark degrees of freedom. Realizing this, we can 
write the complete basis of these effective operators, which can be done most conveniently 
in terms of chiral quark fields, 

{f\nMP\t) = Gj2c,{fi) (/|g,K)(/i) , (15) 

i=l 

where the prefactor G has the dimension of inverse-squared mass, the Cj are dimensionless 
Wilson coefficients,^ and the Qi are the effective operators: 

Q2 = {uLl^,CL) {urYcr) , = {uRit^CR) {urYcr) , 

(16) 

= (ulCr) (urCl) , Q7 = (ulCr) (ulCr) , 

Q4 = (urCl) (urCl) , Qs = (uLCTf^uCR) (uLcr'^'cR) . 

In total, there are eight possible operator structures that exhaust the list of possible inde- 
pendent contributions to |AC| = 2 transitions. Since these operators are generated at the 
scale M where the New Physics is integrated out, a non-trivial operator mixing can occur 
when we take into account renormalization group running of these operators between the 
scales M and /i, with being the scale where the hadronic matrix elements are computed. 
We shall work at the renormalization scale fi = rric ^1.3 GeV. This evolution is determined 
by solving the RG equations obeyed by the Wilson coefficients, 

. (17) 

where 7 represents the matrix of anomalous dimensions of the operators in Eq. f|T6|) (note 
the transposition). Eq. (fT7|) can be solved by transforming to the basis where the transpose 



Throughout this paper, we shaU denote Wilson coeiEcients for |AC| — 1 operators as {c^} and those for 
|AC| = 2 operators as {Ci}. 
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of the anomalous dimension matrix is diagonal, integrating, and then transferring back to 
the original basis Cj. At leading order, we have 



C{fi) = U{fi,M)C{M) , 
where U{n, M) is the evolution matrix, obtained from Eq. f|T7|) by 



f/(/.i,/.2) = V^ r(/.i,/i2r^'''V2/3o v-1 



D 



(18) 



(19) 



In the above, 7''°^ is the vector containing the diagonal elements of the diagonalized trans- 
posed matrix of the anomalous dimensions 7^, is the matrix that diagonalizes 7-^ and 



(20) 



as(/i2) 

For completeness, we display the matrix of anomalous dimensions at leading-order (LO) in 
QCD [151, 



/ 



7 



6 



\ 

















Nc 




12 



-6iV, + ^ 

























6-6iV, + |- 
-24 - ^ 









Nc 









6 + 2iV,-^ 












6- ^ 

















U UiVc jy^ 2 Nc 

_24 — « I 2 TV 



We note that Ref. [15|] also includes the next-to-leading order (NLO) expressions for the 
elements in the anomalous dimensions matrix. However, we perform our calculations at LO 
here since the NLO corrections to the matching conditions in the various models of New 
Physics have generally not been computed. 

Due to the relatively simple structure of 7, one can easily write the evolution of each 
Wilson coefficient in Eq. f[T^ from the New Physics scale M down to the hadronic scale /i, 
taking into account quark thresholds. Corresponding to each of the eight operators {Qi} 
(i = 1, . . . , 8) is an RG factor rj(/i, M). The first of these, ri(/i, M), is given exphcitly by 



^ ' \as{mt) J \as{mb) J 



6/23 



. a;,(/i) ^ 



6/25 



(21) 



M(TeV) 


ri{mc,M) 


r2[mc,M) 


/ 71 ^\ 


/ 71 J'\ 

r4{mc,M) 


/ 71 /r\ 


1 


0.72 


0.85 


3.7 


0.41 


2.2 


2 


0.71 


0.84 


4.0 


0.39 


2.3 



TABLE I: Dependence of the RG factors on the heavy mass scale M. 



and the rest can be expressed in terms of ri (/i, M) as 



r2(/x, M) = 


[ri(/^,M)]i/2, 






r3(/i,M) = 


[n(/i,M)]-\ 


r6(/x,M) 


= n(/i,M) 


r4(/i, M) = 


[ri(/i,M)](i+v^)/6 , 


r7(/i,M) 


= r4(/i,M) 


r5(/x,M) = 




r8(/i,M) 


= r5(/x,M) 



The RG factors are generally only weakly dependent on the NP scale M since it is taken 
to be larger than the top quark mass, mj, and the evolution of as is slow at these high 
mass scales. In Tabled we display numerical values for the ri{fi,M) with M = 1,2 TeV 
and /i = rric — 1.3 GeV. Here, we compute using the one-loop evolution and matching 
expressions for perturbative consistency with the RG evolution of the effective hamiltonian. 



B. Operator Matrix Elements 

We will need to evaluate the D^-to-D^ matrix elements of the eight dimension-six basis 
operators. In general, this implies eight non-perturbative parameters that would have to 
be evaluated by means of QCD sum rules or on the lattice. We choose those parameters 
(denoted by {Bi}) as follows. 



{Qi) = 




(Q5) 


= flMlB, , 


{Q2) = 


-IMMlB, , 


{Qe) 




{Q3) = 




{Q7) 


= -yiMiB 


(Q4) = 




{Qs) 





(23) 
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where (Qj) = {D'^\Qi\D'^) , and fo represents the D meson decay constant. By and large, the 
compensatory 5-factors {Bi} are unknown, except in vacuum saturation and in the heavy 
quark hmit; there, one has Bi ^ 1. 

Since most of the matrix elements in Eq. (1251) are not known, we will need something more 
manageable in order to obtain numerical results. The usual approach to computing matrix 
elements is to employ the vacuum saturation approximation. However, because some of the 
i^-parameters are known, we would like to introduce a 'modified vacuum saturation' (MVS), 
where all matrix elements in Eq. (l23l) are written in terms of (known) matrix elements of 
{V - A)x{V - A) and {S - P) x {S + P) matrix elements So and b'^\ 



(24) 



where we denote A^^c = 3 as the number of colors and, as in Ref. |9|, define 

SE? - . . (25) 

as well as 

(26) 

-DD 

In our numerical work, we take Bf^ = Bd = 0.82, which is the most recent result from 



attice calculation 



IGj, and use the CLEO-c determination Jd = 222.6 ± 



the quenched 
lQ.7+li MeV 

parameters defined in Eq. (125]) for the full operator set relevant to D meson mixing. 



171]. We urge the lattice community to perform an evaluation of the {Bi} 



III. STANDARD MODEL ANALYSIS 



Theoretical predictions of xd and ?/d within the Standard Model span several orders of 
magnitude. Roughly, there are two approaches, neither of which give very reliable results 
because rric is in some sense intermediate between the heavy-quark and light-quark limits. 
Consider, for example, AFd as given in Eq. ( JT^ 

aed = ^ m I J d'xT{n^^^\=\x) ^L^^i=^(o)}|D°) . 
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q 




FIG. 1: Loop diagram for D D . 

To utilize this relation, one inserts intermediate states between the |AC| = 1 weak hamilto- 
nian densities 7^|^^l=^ This can be done using either quark or hadron degrees of freedom. 
Let us consider each of these possibilities in turn. 



A. Quark-level Analysis 

The "inclusive" (or quark-level) approach is based on the operator product expansion 
(OPE). In the rric ^ A limit, where A is a scale characteristic of the strong interactions, 
AM and AF can be expanded in terms of matrix elements of local operators [3] of 
increasing dimensions suppressed by powers of inverse charm quark mass. An instructive 
example concerns a recent analysis of the leading dimension D = Q case ^ in which the 
width difference y-o is calculated in terms of quarks (c/Fig. [1]) and the mass difference xd 
is then found from dispersion relations. The calculation is carried out as an expansion in 
QCD, including leading-order 0{a^) and next-to-leading order 0{al) contributions with 

xd = + and yr, = yj^°^ + y^^^^ ■ (27) 

Here, LO and NLO denote only the corrections at that order and not the full quantity 
computed to that order. Because ruh > Md, AFd experiences no 6-quark contribution.^ 
This leaves only ss, dd and sd + d's intermediate states contributing to the mixing diagram 
of Fig. [U Taking = 0, the mixing loop-functions will depend on 2; = m^/m^ ~ 0.006. 
Table HI] examines in detail the loop-functions for AFq and shows the results of carrying out 
an expansion in powers of z. We see that the contributions of the individual intermediate 
states in the mixing diagram are not intrinsically small - in fact, they begin to contribute at 
0{z^). However, flavor cancellations remove all contributions through 0{z^) for AF^, so the 

^ We ignore here the 6-quark contribution to AA/d; its numerical contribution is subleading (IKidV'crfl ~ 
\VusVcs\ ^ 0.22 whereas ^ 1-8 • IQ-^). 
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Int. State 0{z^) 




0{z^) 


ss 


1/2 


-3z 


3z2 


dd 


1/2 








sd + ds 


-1 


3z 




Total 












TABLE II: Flavor cancellations in AFd- 



net result is 0{z^). Charm mixing clearly experiences a remarkable GIM suppression! The 
corresponding result for AMd turns out to be O^z"^). Summarizing, the leading dependences 
in z for the dimension six contributions are 



(LO) 3 
1/D ~ ^ 



(LO) 2 



(NLO) 2 



(NLO) 2 



(28) 



The source of this ^-dependence is understood as follows. The mixing amplitude is known 
to vanish in the = = limit, so the breaking of chiral symmetry and of SU(3) 
flavor symmetry play crucial roles. Thus, a factor of comes from an SU{3) violating 
mass insertion on each internal quark line and another from an additional mass insertion 
on each line to compensate the chirality flip from the first insertion. This mechanism of 
chiral suppression accounts for the z'^ dependence of x^'~'\ y^^^'^'^ and The case of 

y^'^^ requires yet another factor of to lift the helicity suppression for the decay of a 
pseudoscalar meson into a massless fermion pair. 

Let us next display the LO expressions for y^ and xd (to leading order also in z) j^. 



„(LO)[23] 



y2 ^2 £2 



D 



2 — 2ciC2 — 3Ci 

5 



(S) 



clB^ 



4 



(4 - 2ciC2 - 3ci)B\^ 



(29) 



where C,s = KtsKl '^1,2 are the relevant Wilson coefficients. For our numerical computa- 
tions, we adopt the values used in Ref. [9|], 



rUc = 1.3 GeV 



ci = -0.411 , C2 = 1.208 



a. 



0.406 



(30) 



Numerical results for the LO and NLO contributions, where a discussion of the NLO effects 
can be found in Ref. j^, (c/ Table IIIip reveal that y^ is given by j/nlo to a reasonable 
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LO 


NLO 


LU + iNLU (Central Values j 


2/D 


-(5.7- 


9.5) 


10-8 


(3.9 - 


9.1) 


10-7 


~ 6 • lO-'^ 




-(1.4- 


2.4) 


10-6 


(1.7- 


3.0) 


10-6 


~ 6 • lO-'^ 



TABLE III: Results at dimension-six in the OPE. 



approximation (due to the z dependence discussed above) whereas xd is greatly affected by 
destructive interference between xlo and xnlo- The net effect is to render ?/d and xd of 
similar small magnitudes, at least through this order of analysis. 

The quark- level prediction of xd and ?/d just described is a result of expanding in terms 
of three 'small' quantities, z, A/rric, and as- As a consequence, the use of an OPE to 
describe charm mixing is not entirely straightforward because terms suppressed by higher 
powers of rric could nevertheless be important if they contained relatively fewer powers of m<,. 
However, at the next orders in the OPE one encounters 0{z^^'^) corrections multiplied by 
about a dozen matrix elements of dimension-nine operators and 0{z) corrections with more 
than twenty matrix elements of dimension-twelve operators. This introduces a multitude 
of unknown parameters for matrix elements that cannot be computed at this time. Simple 
dimensional analysis [lol suggests the magnitudes xd ~ ~ 10"^, although order-of- 
magnitude cancellations or enhancements are possible.^ 



B. Hadron-level Analysis 

The D meson mass is not very large, so one might question whether the OPE approach 
discussed in the previous subsection can successfully describe D^-D^ mixing. This is espe- 
cially so since the leading contribution in the SU (3)-breaking parameter enters only as 
a A^/ suppressed contribution in the 1 / m-c expansion, which implies that one has to deal 
with a large number of unknown operator matrix elements. 

As an alternative, one might consider saturating the correlation functions of Eq. f|T4|) with 

^ Any effect of liigher orders in l/nic or asirric) which could produce a z" contribution in the lowest possible 
power n — 1 could yield a dominant contribution to the prediction of and j/d i2o[ [2l] | . Although the 
BaBar and Belle observations of y ^ IQ-^ could be ascribed to a breakdown of the OPE or of duality, it 
is clear that such a large value of ?/d is by no means a generic prediction of OPE analyses. 
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exclusive hadronic states, switching to a purely hadronic description. This approach should 
be valid as the mass of the D meson lies in the middle of a region populated by excited 
light-quark states. In principle, this "exclusive" (or hadronic) approach should sum over all 
possible intermediate hadronic multiplets. Since one has to deal with off-shell hadronic states 
in the calculation of xd, some modeling is necessarily involved. By contrast, a calculation 
of ?/D in this approach is less mo del- dependent. The usual approach to computing xd is to 
first calculate j/d and then use a dispersion relation to obtain xt>. This is appropriate, as 
the contribution due to 6- flavored intermediate states (which appears in xd but not ?/d) is 
negligibly small. 

One possible approach would be to select a set of, say, two-body intermediate states ^, 
and write their contribution to mixing in terms of charged pseudoscalar (P+P^) branching 



fractions 



23 



Vd ^ ^ = ^[DO-^K+K-] + B[D"^TT+TT-] — cos 6ktt BiDO~^K-tt+] ' ^IDO^K+n-] , (31) 

where Sxn is as in Eq. f|T2|) . One can use available experimental data on two-body branching 
ratios to estimate their contribution to ?/d- A dispersion relation then relates ?/d to xd- How- 
ever, the example above explicitly shows the cancellations between states that are present 
within a given SU{3) multiplet. Such cancellations make this procedure very sensitive to 
experimental uncertainties. One would need to know the contribution of each decay mode 

y not feasible at this time. Another possi- 



with extremely high precision, and that is simp 
bility is to model |AC| = 1 decays theoretically 25|]. In this reference, AFd was determined 
in this manner and the result ?/d — 10^'^ was found. This result is, however, smaller than 
the recent BaBar and Belle observations. 

Clearly, is not sufficiently light for its decays to be dominated by just two-body 
final states. Multiparticle intermediate states must also be taken into account in D^-^ 
mixing calculations. In doing so, it is convenient to calculate the contribution of each SU (3) 
multiplet separately, as SU{3) symmetry produces substantial cancellations among members 
of the same multiplet as we saw above. This can be thought of as a long-distance version 
of the GIM mechanism. The surviving contribution is expected to be of second order in the 



^ The simplest intermediate state is a single-particle resonance contribution. Preliminary estimates of 
resonance contributions to D^-D^ mixing appear to be small [2j|, although much remains to be learned 



about the resonance spectrum in the vicinity of the mass. 



15 



S'f/(3)-breaking parameter 20|]. Denoting by ypj^ a value that y would take if elements of 
the final state F belonging to SU (3) representation i?, or Fr, were the only channels open 
for D-decay, one can write y-Q as a sum over all possible Fr's weighted by the D-decay rate 
to each representation, 



Vd 



n) 



(32) 



It is possible to show that yp^ can be computed as [20|] 

^^"^ i:n^F^{D^\nM{n\n^\D^) ' ^ ' 

It should be noted that in the limit of CP-conservation and retaining phase space differ- 
ences as the only source of 5'f/(3) breaking [i.e. neglecting 5'f/(3) breaking in the matrix 
elements), can be computed without any hadronic parameters. This is an appropri- 
ate approximation, as the main contribution comes from the multiparticle (four-particle) 
intermediate state multiplets. For those states, there are multi-kaon modes which are kine- 
matically forbidden. In such cases, phase space effects alone can provide enough SU{?i) 



violation to induce yr, ~ 10" 



20|. In other words, such large effects in y-^ appear for decays 



near the D threshold, where an analytic expansion in SU (3) violation is no longer possible. 
It is interesting that such effects from multiparticle states are not reproduced in the OPE 
calculation, as the resulting contribution does not come from short-distances. 

The use of a dispersion relation for xd then suggests it would receive contributions of 



a similar order of magnitude as those for y-^ [26|. An important difference between the 
resulting values of xd and y-Q is that even retaining phase space differences as the sole 
contributor to SU (3) breaking does not insure cancellation of the hadronic matrix elements. 
However, with some reasonable model-dependent assumptions, one arrives at the conclusion 
that xd ~ l/D ~ 1% [26]. It is thus reasonable to believe that the observed D^-D^ mixing is 
reflecting Standard Model contributions. 



C. Comments 



The above discussions show that, contrary to the B system. Standard Model estimates 
of xd and y^ for the charm system contain significant intrinsic uncertainties. On the other 
hand, SM values near those found by BaBar and Belle cannot be ruled out. Therefore, it will 
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be difficult to attribute a clear indication of New Physics to D^-D^ mixing measurements 
alone. This means that the only robust signal of New Physics in the charm system would 
be the observation of large CP violation, which we will not consider here. Nonetheless, a 
thorough analysis of indirect New Physics contributions is of value, and we find that large 
regions of parameter space can be excluded in many models, placing additional restrictions 
on model building. This will be useful in conjunction with corresponding direct searches for 
New Physics at the LHC 

In what follows, we will take the approach that the New Physics contributions cannot 
exceed the la experimental upper bound for xb- Keeping in mind that this upper limit is 
likely to change as data samples increase and analyses mature, we also display the effects 
of xd < (15.0,8.0,5.0 and 3.0) ■ 10~^ on the parameter space of New Physics scenarios. 
These values are to be used as a guide for how our resulting constraints may change in the 
future. In addition, we will neglect the errors on the determinations of the D meson decay 
constant and B-factors; this will have a small effect on our results given the present large 
uncertainty in the experimental determination of D^-T^ mixing. In all cases, we will neglect 
the possibility of interference between the SM and New Physics contributions. We now turn 
to the examination of various scenarios for physics beyond the SM. 

IV. EXTRA FERMIONS 

The quark sector of the Standard Model can be modified in several ways, and new fermions 
are predicted to exist in many extensions of the SM. They can be classified according to their 
electroweak quantum number assignments; here we consider the possibilities of a sequential 
fourth generation quark doublet (Sect. HV^ . heavy quark iso-singlets (Sect. HVB]) and non- 
SM quarks associated with Little Higgs models (Sect. IIV Cp . The contributions of such 
heavy quarks can remove the efficient GIM cancellation inherent in the short distance SM 
computation and can give rise to D^-D^ mixing at the level of the current experimental 
limit. 
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A. Fourth Generation Quark Doublet 



A simple extension to the Standard Model is the addition of a fourth family of fermions. 
Precision electroweak data severely constrains this possibility. The Particle Data Group [isl 
quotes a restriction on the number of families to be Np = 2.81 ± 0.24 from the oblique .S* 
parameter 



27( 1 alone. We note, however, that the LEP Electroweak Working Group [2^ 



allows for a more generous range of the S parameter from their electroweak fit. In either 
case, this restriction can be relaxed by allowing the T parameter to vary as well, or by 
adding other sources of new phy sics which would participate in the electroweak fit such 
as an extended Higgs sector [29!]. The requirement of anomaly cancellation implies the 
existence of a fourth lepton family as well (almost degenerate to satisfy the Ap constraint 
with m^^ > Mz/2) or an extra right-handed quark doublet. Direct collider searches by the 
CDF and DO Collaborations at the Tevatron currently place a bound 13| on the mass of a 
charged —1/3 fourth generation quark b' of m^/ > 128 , 190 , 199 GeV if the 6'-quark decays 
respectively via charged current interactions into leptons + jets, via FCNC with b' —>■ bZ, 
or is quasi-stable. We recall that perturbative unitarity considerations |30|] in FF — > FF 
scattering restricts the mass of sequential heavy fiavors to be m^? ^ 500 GeV. 

Here, we review the contribution of a fourth generation of quarks to D mixing, keeping 
in mind that some other New Physics may also be present in order to evade the precision 
electroweak constraints and that it also may or may not contribute to the mixing. The 
primary motivation for this discussion is to set up the formalism that will be used in the 
following Sections. 

The Q = —1/3 fourth generation quark contributes to D mixing via a box diagram which 
also contains the SM bosons. Note that since the 6'-quark is not kinematically accessible 
in charm-quark decay, it will not contribute to the dispersive amplitude for xd in Eq. (fT^ . 
The I AC I = 2 hamiltonian at the b' mass scale in the fourth generation model is 31] 

'Hith= ^ ^'^ ^ AjAj5'(xi,Xj)Qi , (34) 



where S{xi the well-known Inami-Lim functions [32[ (given in the Appendix), x. 



mg./MwY, Aj = V*^Vui, and the sum runs over the internal quark fiavors. As discussed in 
the previous Section, there is a strong GIM cancellation in D meson mixing, which leaves 
a sizable contribution only from the heavy b' quark and sets i = j = b' in the above sum. 
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Performing the RG evolution, we obtain at the scale rric 

GlMl 



47r2 



Afe, S {xb' , Xb')ri {rric, Mw)Qi 



(35) 



which in turn gives 



G'pM^ff,MD 



BDXl,ri{mc,Mw)S{xb',Xb' 



(36) 



It should be noted that for ri(mc, Mw), only contributions below Mvi/ are required. 

The value of function of the CKM mixing elements is displayed in Fig. [2] for 

various values of the 6'-quark mass. We see that the la experimental limit of a;D < 11.7x 10^^ 
places sizable constraints in the 6'-quark mixing-mass parameter space. We also show the 
exclusion contours for possible future experimental bounds of xd < (15.0 , 8.0 , 5.0 , 3.0) x 10^^ 
(corresponding to the blue dashed, red dashed, cyan dotted, and green dot-dashed curves, 
respectively) as discussed in the Introduction. We note that the present constraints on the 
CKM mixing parameters |V^fe'V^fe/| ^ 0.002 are an order of magnitude stronger than those 



obtained from unitarity considerations 



13| of the CKM matrix. 




, I , , , , I , 

0.002 0.003 
|V„b-Vcb'l 




0.002 
IVnb'Vcb'l 



FIG. 2: Left: xd in the four generation model as a function of the CKM mixing factor |V^{,/Kj6'| for 
6'-quark masses of 200, 300, 400, and 500 GeV from bottom to top. The la experimental bounds 
are as indicated, with the yellow shaded area depicting the region that is excluded. 
Right: The present la excluded region in the mass-mixing parameter plane, as well as possible 
future contours taking xd < (15.0,8.0,5.0,3.0) x 10"'^, corresponding to the blue dashed, red 
dashed, cyan dotted, and green dot-dashed curves, respectively. 
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B. Heavy Vector-like Quarks 



The next possibility of interest is the presence of heavy quarks which are SU{2) l singlets 
(so-called vector-like quarks) 33|. We will consider both charge assignments Q = +2/3 
and Q = —1/3 for the heavy quarks. Both choices are well motivated, as such fermions 
appear explicitly in several models of physics beyond the Standard Model. For example, 



351], with one for each of the 



weak isosinglets with Q = —1/3 appear in Eq GUTs 
three generations {D, S, and B). Weak isosinglets with Q = +2/3 occur in Little Higgs 



theories [36|, |37| in which the Standard Model Higgs boson is a pseudo-Goldstone boson, 
and the heavy iso-singlet T quark cancels the quadratic divergences generated by the top 
quark in the mass of the Higgs boson. 



1. Q = —1/3 Singlet Quarks 

We first consider the class of models with Q = —1/3 down-type singlet quarks. For this 
case, the down quark mass matrix is a 4 x 4 array if there is just one heavy singlet (or 6 x 6 
for three heavy singlets as in Eq models). As a consequence, the standard 3x3 CKM matrix 
is no longer unitary. Moreover, the weak charged current will now contain terms that couple 
up-quarks to the heavy singlet quarks. For three heavy singlets, we have 

L^t^ = ^V,^W''u,,Ll>.D^ , (37) 

where Ui^L = {u, c, t) l and = {D, S, B) refer to the standard up quark and heavy isosinglet 
down quark sectors. The {Via} are elements of a 3 x 6 matrix, which is the product of the 3x3 
and 6x6 unitary matrices that diagonalize the Q = +2/3 and Q = —1/3 quark sectors, 
respectively. The resulting box diagram contribution to AMd from these new quarks is 
displayed in Fig. [31 Assuming that the contribution of one of the heavy quarks (say the S 
quark) dominates, one can write an expression (similar to that in Eq. (1361) ) for xd |31[ |. 

BDiy^s^us) ri{m^,Mw)f{xs) , (38) 

where xs = {ms/My/Y and f{xs) xs{l + 61n(a;5)) for large xs- The light-heavy mixing 
angles jV^^^V^sl^ should go as l/m^ for large mg to keep the contribution under control. 



The current bound on 



2 

Ks^sl from unitarity of the CKM matrix is not very stringent. 



iKs^sl < 4 X 10 ^ 13|. An S'-quark mass in the range 0.2 to 1 TeV gives rise to a mixing 
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FIG. 3: Box contribution from heavy weak-isosinglet quarks. 

contribution that can exceed the current experimental hmit in Eq. (jlj). Hence a singlet 
heavy quark of charge —1/3 can give rise to Xd near the current experimental limit. 



In the -Ee-based model proposed by Bjorken et al [38|], the 6x6 mass matrix has an 
especially simple form. The resulting 6x6 mass matrix has a pseudo-orthogonality property 
which implies that the 3x3 CKM matrix, although not unitary, satisfies 

E (^ckm):, {VcKuh, = . (39) 

i=l 

The analog of this condition in the up quark sector does not hold, and as a result, there are 
no new FCNC effects in the down quark sector. For D^-D^ mixing, the prediction is now 
that (recall capital lettering is used to denote the heavy quark) 

KsVuS? = sl\V:,Vus\' ^ sl\^ , (40) 

where |V],*Kis| — A ~ 0.22 and S2 is the (small) mixing parameter describing the mixing 
between the light s quark and the heavy S quark. Using the experimental values in Eq. (jlj), 
we can place bounds on S2 for a given mass ms-, e.g. S2 — 0.0009 for rris — 0.5 TeV. We 
present the constraints on S2 vs. ms in Fig. HI 



2. Q = +2/3 Singlet Quarks 

Next, consider the possibility of weak isosinglet quarks having charge Q = +2/3. These 

are present in some theories beyond the SM, including for example. Little Higgs models 

which will be discussed in more detail in the following section. Here, we present the general 

formalism for this scenario. 

The presence of such quarks violate the Glashow-Weinberg-Paschos naturalness condi- 

□ 

tions for neutral currents [391 . Since their electroweak quantum number assignments are 
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FIG. 4: Left: xd in the singlet Q = —l/2> quark model as a function of the singlet quark mass for 
various values of the mixing angle S2 = 0.0001 , 0.0002 , 0.0003 corresponding to the solid, dotted, 
and dashed curves, respectively. The la experimental bounds are as indicated, with the yellow 
shaded area depicting the region that is excluded. 

Right: The present la excluded region (short-dashed curve) in the mass ms and mixing angle S2 
parameter plane for the singlet Q = —1/3 quark in the model of Bjorken et al. [38] described in 
the text. Possible future contours are also shown, taking xd < (15.0,8.0,5.0,3.0) x 10^'^ from 
top to bottom, corresponding to the solid, medium dashed, long dashed, and longer dashed curves, 
respectively. 

different than those for the SM fermions, flavor changing neutral current interactions are 
generated in the left-handed up-quark sector. Thus, in addition to the charged current 
interaction 

= A^V^^^^liA^lW^' . (41) 



there are also FCNC couplings with the boson 331], 



& = ,J , A.,«,„^7.«..L^°^ . (42) 

Here, Vai is a 4 x 3 mixing matrix with a running over 1 ^ 4, i = 1 — > 3, and with the 
CKM matrix comprising the first 3x3 block. In this tree-level contribution to 

AMe) is generated from Z°-exchange as shown in Fig. [5l This is represented by an effective 
hamiltonian of the form 

^^/^ = Scos'g.Ml ^-7.c.^.7^c, , (43) 

where unitarity demands 
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FIG. 5: Tree-level contribution from Z'^-exchange. 



Taking the la ranges for the experimentally determined values of the CKM elements 13 1 
yields the constraint Xuc < 0.02. This hamiltonian is just a particular case of a more general 
relation (Eq. f H9|) ) appearing in the following section.^ The QCD running from n = Mz to 
jj, = rric for this effective hamiltonian is trivial, leading to 

-,2 



9_ 

where it should be noted that for ri{mc., Mz)-, only contributions below Mz are required. 
This hamiltonian leads to 

(2/3) "^GpfliM-D j3 x2 / N 

= ^^Yd n{m^,Mz) . (46) 

We present this contribution to xd from models with a singlet Q = 2/3 quark in Fig. [61 Note 
that the bound on the mixing X^c from D^-T^ mixing is roughly two orders of magnitude 
better than that from unitarity constraints of the CKM matrix. 

C. Little Higgs Models 



Little Higgs models [36|, |37|, |40|] feature the Higgs as a pseudo Nambu-Goldstone boson 
of an approximate global symmetry that is broken by a vacuum expectation value (vev 
hereafter) at a scale of a few TeV. This approximate symmetry protects the Higgs vev 
through one-loop order relative to the ultra-violet (UV) cut-off of the theory which appears 
at a higher scale. The breaking of this symmetry is realized in such a way that the Higgs 
mass only receives quantum corrections at two-loops. In these models the one-loop quadratic 
divergent contributions to the Higgs mass in the SM are canceled by a new particle of the 
same spin. These models thus predict the existence of new charged Q = +2/3 vector- like 
quarks, gauge bosons, and scalars at the TeV scale. 



The specific correspondence is Cn — 0,Cl = Cy^ — g^A^^/(4cos^ ^m) and Mz' — Mz- 
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FIG. 6: Value of xd as a function of the mixing parameter A^c in units of 10-2 in the Q = +2/3 
quark singlet model. The la experimental bounds are as indicated, with the yellow shaded area 
depicting the region that is excluded. 

The most economical model of this type, in that it introduces the minimal number of new 



fields, is known as the Littlest Higgs [36|, l37| . It is based on a nonlinear sigma model with 
SU{5) global symmetry that is broken to the subgroup 5*0(5) by a vev /. / is generated 
by strongly coupled physics at the UV scale A ~ iirf ~ 10 TeV. The fourteen Goldstone 
bosons remaining after this symmetry breaking yield a physical doublet and complex triplet 
under SU{2), which remain massless at this stage. The SU{5) contains a gauged subgroup 
[SU (2) X U (1)]2 which is also broken by / to the SM electroweak gauge group. The remaining 
four Goldstone bosons are then eaten by a Higgs-like mechanism and give mass, of order 
/ ~ 1 TeV, to the gauge fields of the broken subgroups. Masses for the complex triplet are 
generated at the TeV scale by one-loop gauge interactions. The neutral component of the 
complex doublet plays the role of the SM Higgs, which receives its mass at two-loops from a 
Coleman- Weinberg potential, giving /i^ ~ /^/IGtt^. Thus the natural scale for / is around 
a TeV; if / is much higher, the Higgs mass must again be finely tuned and this model no 
longer addresses the hierarchy problem. 

The minimum physical spectrum of this model below a TeV is thus that of the SM with 
a single light Higgs. At the TeV scale, there are four new gauge bosons (an electroweak 
triplet and singlet), the scalar triplet, and a single Q = +2/3 vector-like quark T. Other 
variants of Little Higgs models may expand this particle content at the TeV scale. 
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In general, the new vector-like T quark can contribute to D^-D^ mixing. Since it has 
different electroweak quantum numbers from the Q = +2/3 quarks in the SM, FCNC in- 
teractions will be induced in the left-handed up quark sector. This generates a tree-level 
Z boson exchange contribution to D mixing as depicted in Fig. [5l This was considered in 
Ref. where a specific ansatz for the 4 x 4 up quark mass matrix was employed, leading 
to a quite small contribution to AM^. In general, however, one expects the quark mixing 
to be of order v/f and the contribution of the flavor-changing Z interaction induced by the 
existence of the T quark can be sizable, as discussed in the previous section. Current data 
can be used to constrain the mass and mixing of the T quark, and the results in Fig. [6] are 
applicable in this case. Tree-level contributions to D mixing with the exchange of the new 
heavy neutral gauge bosons will likewise be generated, as the fermionic coupling for these 
fields is also proportional to the fermion's third component of weak isospin. The operator 
structure is exactly the same as that given in the previous section, but the magnitude of 
these contributions will be suppressed relative to the SM Z boson exchange by the heavy 
mass of the new neutral gauge bosons. In addition, there could also be new contributions 
to rare D meson decays, as discussed recently in the work of Chen et al. 421]. 



It has been shown 



43| that a global fit to the precision electroweak data set places a 



significant limit, which is roughly parameter independent, on the vev / in the Littlest Higgs 
model of / ^ 4 TeV. Variants of this model, employing different global symmetries, can 
reduce this constraint somewhat 4J| . In addition, a discrete symmetry, called T-parity, can 
be introduced 45|] to alleviate the electroweak bounds. This symmetry is analogous to R- 
parity in supersymmetry and has the consequence that T-parity is conserved in interactions 
and that the lightest T-odd particle is stable. This provides a natural Dark Matter candidate 
in these models. Recent work by Hill and Hill 46|| has shown that anomalies in topological 
interactions can break this discrete symmetry and thus T parity is no longer an exact 
symmetry; the resulting phenomenology has yet to be worked out. However, the Ultra- 
violet (UV) completion of the theory may or may not allow for the terms which break 
T-parity, and thus a general statement on the presence of this discrete symmetry in the 
low-energy theory cannot be made. 

In addition to the tree-level contribution discussed above. Little Higgs models with T 
parity can give rise to a loop contribution to D^-D^ mixing involving the exchange of the 
heavy gauge bosons and new mirror fermions, which are present in this form of the model. 
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In fact, three vector-like doublets of mirror fermions are introduced in Little Higgs models 



45| . The effective hamiltonian 



with T parity in order to evade compositeness constraints 
relevant to D^-D^ mixing for this contribution is 47 1 

n.n = ^^EePer^.H(^...)Qi . (47) 

Here, ^^^^ corresponds to the relevant elements of the weak mixing matrix in the mir- 
ror fermion sector which parameterizes the flavor interactions between the SM and mir- 
ror fermions. The quantity (given in the Appendix) is the loop function computed in 



Ref. [47[]; it depends on Zi = ml^jM^^, where rriMi is the mass of the i^^ mirror quark 
doublet and Wh represents the heavy charged gauge boson mass. The Qi operator appears 
since the heavy gauge bosons Wh have purely left-handed interactions. The RG running of 
this hamiltonian is trivial and leads to a factor of ri{mc, M). The resulting contribution to 
the mass difference is 

— 24^r2f^ ^^jiZ^^i Q 'FH{zi,Zj)ri{m^,M) . (48) 



This has recently been computed in Ref. [48| in light of the recent experimental measurement 
of D^-D^ mixing, where it is found that this contribution can saturate the experimental 
bounds. 

In Little Higgs models with T parity, there is an additional tree-level contribution arising 
from the interaction vertex q T^~^ where Zfj represents either of the heavy neutral gauge 
fields and T^~^ is the odd T-parity quark. T*-"-* couples to the weak eigenstate of the T- 
parity even quark T^'^\ which receives its mass from the same Q = +2/3 Yukawa term that 
is responsible for the up quark masses. This induces mixing between the quarks in the up 
quark sector, resulting in FCNC interactions of the SM quarks with the exchange of the new 
heavy neutral gauge fields. The generic formalism for this contribution will be discussed in 
the next section, however this particular contribution is thought to be small \^\ - 



V. EXTRA GAUGE BOSONS 



Many theories with physics beyond the SM have extended electroweak gauge symmetries, 
whose hallmark are the existence of new heavy neutral and charged gauge bosons. We note 
that scenarios with extended gauge symmetries also generally contain new fermions which 
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are required for anomaly cancellation, as well as an extended Higgs sector to facilitate the 
extended symmetry breaking. The additional heavy gauge bosons are produced directly 
at hadron colliders via the Drell-Yan mechanism, and the search limit on the masses of 
new W, Z' gauge bosons with SM couplings is approaching 1 TeV from Run II data at the 
Tevatron. The lower bound on the mass of a SM-coupled Z'{W') is 923(965) GeV from 
CDF (DO) [49!. The LHC will be able to search for these particles with masses up to 5 
TeV [50] . There are many such models that yield large FCNC effects in the up-quark sector. 



A. Generic Z' Models 



It is possible that a new heavy Z' boson has flavor changing couplings in the up quark 
sector. Here, we examine a generic tree-level FCNC interaction that mediates D mixing via 
Z'-exchange, analogous to the transition depicted in Fig. [5l While the discussion presented 
here is quite general, many string-inspired models have extra U(l) gauge symmetries that 



35, 



lead to extra Z'-bosons with possible flavor- changing couplings 
The effective four-fermion hamiltonian just below the Z' scale is 

1 
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52|. 



TCz' 



2M|, 



Cl Qi + 2ClCr Q2 + Cr Qe 



(49) 



where the dimensionless flavor changing couplings Cl.r are the model dependent inputs to 
the calculation. This is the most general effective hamiltonian and assumes flavor-changing 
interactions occur in both the left- and right-handed sectors. We first perform a general 
analysis and will then consider some particular occurrences of a Z' in the Sections below. 
We introduce the Wilson coefficients Ci^2 fi{Mz') by matching at the Z' mass scale. 



C2(MzO = 2ClCr , C.iMz') = Cr2 



(50) 



with all other Wilson coefficients being zero at this scale. Assuming that Mz> > rrit and 
performing the RG running of Eq. (jUj), we obtain the effective hamiltonian at the scale 
fi = rUc, 

1 



2Ml, 



Ci(mc) Qi + C2(mc) Q2 + C3(mc) + C6(mc) Qt 



(51) 



with 



Ci(m,) =ri(m„MzOCi(MzO , 
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CaK) = r2{m„ Mz')C2{Mz') , 

Cs{m,) = ^ hK, Mz') - rs{m„ Mz>)] C2{Mz') , 

C,{m,)=r,{m,,Mz')C,{Mz') . (52) 

The presence of Q3 in Eq. (ISTl) is due to operator mixing in the RG running. Finally, as 
a check note that for the case of no evolution (r^ 1) we obtain the expected behavior 
Q(/i) ^ Q{Mz'). 

Upon evaluating the D^-to-D^ matrix elements, we obtain the Z' tree contribution to xd, 



(z') _ fpMp 



^ (Ci(m,) + C,{m,)) - C2{m,) + ^) 



(53) 



where we have made use of Eqs. fl24l). fl26l) . 

Eq. (!53l) can be used to relate the input parameters of Z' models (Cl and Cr) to some 
value of xd- Taking xd < 11-7 x 10"'^, particularly simple expressions are obtained for the 
limiting cases: 

1. Cr = (the case with Cl replaced by Cr yields identical limits): 



Mz' fflMnBnnimc, Mz') \ 



Cl \ S^dEd 

2. Cl = Cr = C: 

Mz> (flMoBDKV'^ 



> 8.9 X 10^ GeV , (54) 



> 3.4 X 10^ GeV , (55) 



c V 2xdE,5 ; 

where k= \Ari{m^, Mz')/3 - 8r2im^, Mz')/9 - {1 + Q7])rs{m^, Mz')/9\. 

In the above, we have fixed the slowly varying ri{mc, Mz') and K{mc, Mz>) at the middle 
of their ranges. The results in Eqs. fl53l) - fl55l) which constrain combinations of the flavor- 
changing couplings and the Z' mass (see Fig. [7]), can be applied to fit the needs of the 
NP model builder. For example, in Eq. fl5^ the choice Cl = 1, Cr = implies Mz' > 
8.9 X 10^ TeV or ahernatively taking Mz' = 1 TeV yields'^ the bound Cl < 1.1 x 10"^. 



After the work described in this section was completed, we received Ref. 53| in which similar results were 
discussed. 
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FIG. 7: xj) in a generic Z' model as a function of the Z' boson mass (normalized by the flavor 
changing coupling constant) for the case Cl = Cr = C. The la experimental bounds are as 
indicated, with the yellow shaded area depicting the region that is excluded. 

B. Family Symmetries 

One class of new physics possibilities to be explored in the TeV range is family (horizontal) 
symmetries. The gauge sector of the Standard Model lagrangian, 

i^SM = ^QCD + CsU{2)yiU{l) + ^^H , (56) 

actually exhibits a very large global symmetry viz. SU (3) l x SU{?>)r x U{1) x f/(3) l x ^^(3)/?. 
The hope is then that some subgroup G of this large symmetry is shared by the Higgs sector 
Ch and the gauge symmetry of the full lagrangian becomes SU{?>)c x SU(?>)l x U{1). The 



group G acts on the families horizontally 5J], and , of course, eventually G has to be broken, 
preferably spontaneously. 

The symmetry G can be implemented locally, so there will be flavor- changing interactions 
mediated by new gauge bosons. The symmetry is broken spontaneously, making the gauge 
boson massive with new scalar fields being introduced in addition to the standard Higgs 
field. 



As a prototype, let us consider a very simple model [55|. We consider the group SU{2)g 
acting only on the first two left-handed families (it may be regarded as a subgroup of an 
SU{?>)g, which is broken). Spontaneous breaking of SU{2)g makes the gauge bosons Gi 
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massive. The LH doublets 



and 



V J 



(57) 



transform as Iq = 1/2 under SU{2)g, as do the lepton doublets 







and 






L 







(58) 



and the right-handed fermions arc singlets under SU{2)g- The superscript refers to the fact 
that these are weak eigenstates and not mass eigenstates. The couplings of fermions to the 
family gauge bosons G is given by 



(59) 



where / denotes the couphng strength and r are the generators of SU{2)g- We define the 

mass basis by 



u,. 



U„ 



U, 



(60) 



In the limit of CP-conservation each of the three 2x2 matrices Ud, Uu, and 1/^ is characterized 
by one angle: 6d, 6u, and (where the Cabibbo angle is 9^ = 0^ — Od)- One then finds for 
the couplings in the fermion mass basis: 

^ = f Gi^|sin26'rf {dLliidL - slI^sl) + sin26'„ {uLlt^UL - cl'^^cl) 
sin 291 (eLT/^eL - T^lIixI^l) + cos 26*^ (dLl,j.SL + SLl^di) 

cos 2eu {uLlijiCL + Cl^hUl) + COS 261; {eLl^ijJiL + l^Llti^L) } 
+ G'3^{cOs26'<i (d.L'^^dL - Sl1,j.Sl) + C0S26'„ {uLl^iUL - Cl^i^Cl) 

+ COS 201 (eL^u^eL - JIlIhI^l) - sin 261^ (dLlnSL + SLlndh) 

- sin 26^ {ul'J^cl + CL^i^UL) - sin 261 {cl'^^hl + TIlIh^l) }] ■ (61) 

For simplicity we assume that after symmetry breaking the gauge boson mass matrix is 
diagonal to a good approximation in which case Gi^ are physical eigenstates and any mixing 
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FIG. 8: xj^ as a function of rrii/f for the gauge boson mass ratios mi/m2 = 0.01,0.5,0.7 
corresponding to the sohd, dotted, and dashed curves, respectively. The la region for xd is also 
shown, with the yellow shaded region depicting the excluded region. 

between them is neglected. This lagrangian clearly introduces tree-level FCNC interactions 
and gives a contribution to D^-D^ mixing of 

2 /cos^ sin^2g. 1 \ _ _ ^ 



A simple symmetry breaking pattern (see, e.g.^ [55|]) leads to m\ = ^ m2. Since the 
effective hamiltonian only involves the operator Qi, the RG-running is simple and leads to 
the following structure at the scale, 

2 2 

T-CFsimc) = /Vi(me, M) % ^ ^ ulI^iClUl'^^cl , (63) 

Tfi^Tfi2 

where M is the smaller of the new gauge boson masses nii and m2. This leads to a value 
for x^^^ of 

^,,,^2^|i^/^^^_»|X ^^^^ 



Using the available experimental data on D^-D^ mixing parameters, yields constraints on 
the masses of the family symmetry-mediating gauge bosons. They are presented in Fig. [H] 
for rrii/f. 



C. Left-Right Symmetric Model 



A puzzling feature of the SM is the left-handed nature of the electroweak interactions. 
A long-standing possible remedy, known as the Left-Right Symmetric Model (LRM) 56| . 
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seeks to restore parity at high energies by enlarging the gauge symmetry to SU{2)l x 
SU{2)r X U{1)b-l- This model can be embedded into an 5*0(10) (or Eq) GUT structure 
which then provides a natural mechanism (seesaw) for generating light neutrino masses. 
A supersymmetric version of a left-right symmetric 5*0(10) GUT model yields the correct 
prediction [57| for ^ 6y^{Mz) and as{Mz), while allowing for the masses of the new 

gauge bosons {Zr and W^) associated with the SU{2)ji symmetry to be of order a few 
TeV or less. Light masses for the new gauge bosons can also be obtained in models with a 
horizontal symmetry 58|. Manifest left-right symmetry dictates that the right-handed gauge 
coupling take on the same value as the left-handed SM coupling ql and that the elements of 
the right-handed CKM matrix be equal to their left-handed counterparts. In this case, the 
direct search for new gauge bosons at the Tevatron places the bound \a^^ of > 788 GeV 
on the mass of the charged right-handed gauge boson. 

The Zji has flavor conserving couplings in this model, and thus does not mediate D^-lf' 
mixing. The charged right-handed gauge field, however, can participate in virtual exchange 
in a box diagram, in association with the SM Q = —1/3 quarks, and gives a contribution to 
meson mixing. In fact, the strongest bound on the mass of the Wn, which is Mr > 1.6 TeV 
in the limit of manifest left-right symmetry, is derived from its contribution to K'^-l{^ 
mixing [s^. 

However, there is no compelling theoretical reason to adopt manifest left-right symmetry 
and in more general models the elements of the right-handed analog of the CKM matrix can 
take on any values, while still respecting unitarity. In addition the ratio of gauge couplings 
can vary [60| between Xw/{1 — Xw) = 0.55 < gn/gi = k < 1 — 2. In this case where 
manifest left-right symmetry is dropped, the bounds from Kaon mixing are softened to 
^ 300 GeV 61| and the direct collider searches are significantly weakened j62]. 

The I AO I = 2 hamiltonian at the right-handed mass scale is given by 

-Hlrm = ['^'K^^cf K"ViV(a:r , (3)Q2 + /.^(K^cf )'5(xf )Q6] , (65) 

with xl = mf /Mf, (3 = M^/Mj^, V^'^ denote the left- and right-handed CKM matrix, and 
the quantities S{x) (an Inami-Lim function) and J{x, (3) are given in the Appendix. The first 
term in this hamiltonian corresponds to the exchange of one Wr and one Standard Model W 
boson in the box diagram, while the second term represents the contribution where only the 
Wr participates. Here, we ignore mixing between the left- and right-handed gauge bosons. 
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The identification of tlie matcliing conditions at tlie liigli scale, 6*2,6 (^r), are obvious. Tlie 
RG evolution to tlie charm scale yields the effective hamiltonian at = rric 

T-thRM = ^jp[C2{mc)Q2 + Ci{mc)Qz + CQ{mc)Q&] , (66) 

where operator mixing has induced the dependence on similar to the generic Z' case 
discussed above, with 



5 



C2K) =r2K,Mfi)C2(MR) 
C3(m,) = ^hK, Mr) - r3(me, Mr)]C2{Mr) , (67) 
C6(mc) = rQ^rric, MR)C6{mc, Mr) . 

Evaluating the hadronic matrix elements yields 

^d''"'^ = IItT?^ [-10C2(me) + 7C3(m,) + 8Ce{m,)] (68) 



upon employing the vacuum saturation approximation and taking rj = 1 in Eq. (!26|l . 

It is clear that in the case of manifest left-right symmetry, the combination of the small 
(KiiKfe*) CKM elements and the suppression will result in a very small value for 

_^gjRM)^ However, it is possible that for non- manifest left-right symmetry, where the right- 
handed CKM elements {Vj^V^;*) may take on larger values, that a significant effect may 
be generated. We examine this scenario, taking {V^V£*) to lie in the range 0.001 — 0.5, 
where 0.5 is the maximum value that this quantity can attain while respecting unitarity of 
the right-handed CKM matrix. Our results are shown in Fig. [9] for various values of Mr, 
taking k = Qr/ql = 1- We see, that even for the most extreme values of the parameters, 
the LRM contribution to D^-D^ mixing never reaches the experimentally determined value. 
Since both of the loop functions J{x,P) and S{x) go as ml, the suppression from the small 
internal quark masses dominates this result. The dip in the curves results from interference 
due to operator mixing. 

This exercise shows that a generic scenario with a new heavy charged gauge boson that 
participates in the box diagram for D mixing will not induce sizable contributions to the 
neutral D meson mass difference, unless it is accompanied by new heavy Q = —1/3 quarks. 



The recently proposed Twin-Higgs models 63|], which are based on a SU{2)r x S{2)l x 



U{1)b-l gauge symmetry, contain an extended top quark sector with heavy Q = +2/3 
quarks and would give significant contributions to K ,Bd^s mixing, but not to AMd. 
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FIG. 9: xd in the Left-Right Symmetric model with non-manifest left-right symmetry as a function 
of the right-handed CKM mixing factor \V^V^^*\ for Mr = 300 , 600 , 1000 , 1500 GeV from top to 
bottom. The la experimental bounds are as indicated, with the yellow shaded area depicting the 
region that is excluded. 

D. Alternate Left-Right Model from Eq Theories 

An alternative to the conventional Left-Right Synametric Model discussed above is pos- 
sible in supersymmetric Eq Grand Unified Theories [35]. This model is also based on the 
low-energy gauge group SU{2)l x SU{2)r x f/(l), but makes use of ambiguous fermion 
assignments within the fundamental representation of Eq [64]. The additional right-handed 
charged and neutral gauge fields in this model have different properties than in the tradi- 
tional scenario. A single generation in Eq theories contains 27, 2-component fermions (in 
contrast to the 16 fermions per generation in SO(IO)), and quantum number ambiguities 
arise that allow the T^l(r) assignments to differ from their customary values for the i'l,r , ^l, 
and dji fields. In the quark sector, the right-handed up-quarks then form SU{2)ji doublets 
with the exotic Q = —1/3 vector singlet quark, Dr, that is present in the 27 representation 
of Eq, while the SU{2)l doublet {u,d) remains unchanged. Di and are then singlet 
fields under all gauge symmetries. This allows, for example, the right-handed W boson 
to couple the right-handed up-quark sector to the singlet quark Dr. Examination of the 
superpotential for this model shows that the Dr takes on the quantum number assignment 
of a leptoquark, while the Wr carries negative R-parity and non-zero lepton number, and 
thus cannot mix with the Wl of the SM or couple to the down-quark sector. The usual 
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constraints on right-handed gauge bosons from the — Ks mass difference and polarized 

/i decay are thus evaded in this scenario. ^ ^ 

These exotic particles can induce significant contributions to D^-D^ mixing [65] via Wr 
and Dfi exchange in a standard box diagram. Note that since the heavy D^i quarks are not 
kinematically accessible in charm-quark decay, there is no dispersive amplitude in this case. 
The interactions of the right-handed W boson take the form 

c = ^v,fua,{i + i,)D,w^R ' ^^^^ 

where i,j are generational indices and V^^ is the right-handed analog of the CKM quark mix- 
ing matrix governing the right-handed charged currents. The effective | AC| = 2 hamiltonian 
at the scale of the right-handed interactions is 

'^ALRM = ^J}m^ E(K^Kf )(KiV;f )^(^. ^,)Q6 , (70) 

where the sum over i,j extends over the three generations of quarks, S{xi the 
standard Inami-Lim functions [32] (given in the Appendix) and Xi = mj^^ jM^j with Mr 
being the mass of the new right-handed gauge boson. Note that this expression mirrors that 
in Eq. except for the presence of the right-handed operator Qq. Matching at the scale 
Mji yields 

C.iMn) = A E (Kf )(K5Kf )^(^., X,) . (71) 
Performing the RG evolution we obtain at the charm quark scale 



nLR = ^C,{m,)Qe , (72) 

R 

with CQ{mc) = rQ{mc, MFi)CQ{Mii). This yields the contribution to the D meson mass 
difference 

^(ALRM) ^ |^^r6(m„ Mn) E(Kf e)(K5Kf )^(^. ■ (73) 

D R ij 

The magnitude of these contributions is determined by the form of the right-handed 
quark mixing matrix, the degeneracy of the 3 generations of quarks, as well as the right- 
handed mass scale. If the quarks are fully degenerate, then a right-handed GIM mechanism 
is operative due to the unitarity of and this contribution to D mixing vanishes. If there 
are mass splittings between the three generations of -Di^^j, then the observed value of D 
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mixing can place bounds on the size of these sphttings. Here, we will examine the case 
where takes on the form of the left-handed quark mixing matrix, i.e., it displays the 
hierarchal structure of the CKM matrix, and derive constraints on this mass splitting as a 
function of the right-handed mass scale. Our results are presented in Fig. (TDl The value of 
function of M/j, the mass of the right-handed charged gauge boson is displayed 
for various mass splittings, Am/m/)^, where Am = — rrij is taken to be constant between 
the first and second as well as second and third generations. We also display the constraints 
the present la experimental bound of xd < 11.7- 10^^, as well as contours for future possible 
values of xd, places on the - mass splitting parameter plane. The results are shown for 
two representative values of the first generation D-quark mass, = 500 and 2000 TeV. 
We see that the mass of the Wr is restricted to be Mr > 1 — 5 = 2.0 for large values of the 
quark mass splittings. 



E. Vector Leptoquark Bosons 

In most recent papers on the subject, scalar leptoquarks (particles with both quark and 
lepton quantum numbers) have usually been associated with i?-parity- violating SUSY scalars 
(to be considered in Sect. IVIII C|) . However, vector leptoquarks (VLQs) are also a possi- 



bility 



661]. For example, they natural 



ly arise in Grand Unified Theories, where quarks and 



leptons belong to the same multiplet 67|] . Many New Physics models where leptoquarks are 
introduced as fundamental vector particles imply that they serve as gauge particles medi- 
ating quark and lepton-number-changing interactions resulting from GUT-model symmetry 
groups. Those symmetries are usually broken at a rather high scale, of order 10^^ GeV, 
which implies that, barring fine-tuning, vector leptoquarks receive masses near the GUT 
symmetry-breaking scale. Yet, some models exist where leptoquarks receive masses at a 
lower scale. In addition, more exotic constructions, such as preon (composite) models, could 
also contain vector leptoquarks. In those models, however, leptoquarks are composite par- 
ticles with masses that are of the order of compositeness scale. Thus, observations of effects 
of VLQs could potentially probe physics at a very high mass scale. It is for these reasons 
vector leptoquarks are searched for experimentally. Collider searches at the Tevatron for 
the direct production of vector leptoquark pairs yield the constraint myiQ > 290 GeV ^] 
from Run I data for second generation leptoquarks which decay into muons. 
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FIG. 10: Left: xd in the alternate left-right model as a function of the mass of the for various 
values of the singlet quark mass sphttings, A/m = 0.1,0.2,0.5,1.0 from bottom to top. The la 
experimental bounds are as indicated, with the yellow shaded area depicting the region that is 
excluded. 

Right: The present la excluded region in the Mji - mass splitting parameter plane, as well as 
possible future contours taking xd < (15.0 , 8.0 , 5.0 , 3.0) x 10""'^, corresponding to the blue dashed, 
red dashed, cyan dotted, and green dot-dashed curves, respectively. The mass of the first generation 
singlet quark Di is taken to be 500 and 2000 GeV as labeled. 

While there are a number of phenomenological studies of vector leptoquarks 0, [t^ , a 
general problem exists with placing constraints on VLQs from indirect measurements, and 
in particular from D^-lf mixing. This is because their couplings are model-dependent. In 
particular, loop calculations with massive composites, i.e. non-gauge leptoquarks, receive 
contributions that are divergent and must be regulated by the compositeness scale. For 
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the case of gauge leptoquarks, and in the absence of a GIM-hke mechanism in the lepto- 
quark box diagram, one can choose a gauge (such as the Feynman gauge) to unambiguously 
compute the effects of leptoquark interactions. In that gauge, however, one also must add 
contributions from unphysical states responsible for the generation of the VLQ masses. In 
a specific model, the interactions of the unphysical states are fixed and their contributions 
are readily computable. However, this then becomes a rather model-dependent procedure 
because VLQ masses can be generated by various means, including some version of the Higgs 
mechanism, or a Frogatt-Nielsen mechanism, etc. Rather than rely on a specific model, in 
what follows, for generality, we shall follow the approach of Ref. [691] and obtain bounds on 
the couplings of gauge VLQs by dropping the contributions from the unphysical states. 

In general, a VLQ could couple to both left-handed and right-handed fermions, so we 
shall assume the general form of the coupling. We note, however, that there are stringent 
bounds [7l| from low-energy data if leptoquarks couple to both left- and right-handed states. 



and it is generally assumed that their couplings are chiral. For a quark of flavor q and a 



lepton £, we adopt the interaction vertex 
to D meson mixing, 

1 



which leads to the contribution 



X 



(VLQ) 
D 



{\ll{Qi)+2\lr{Q2) + \rr{Q&)) 



10 

(All(Q7) + 2Ai^(Q3) + Xrr{Qa)) 

9 m 



LQ 

i^LL + ^Rr) 



(74) 



where Xpp> = (Ap'^Ap"^ (^ApfAp/"^, and we neglect 0{mc/mLQ) corrections in the last 
line. The resulting bounds on VLQ interactions are displayed in Fig. [TTl 



VI. EXTRA SCALARS 



No known physical principle restricts the number of Higgs multiplets that can participate 
in electroweak symmetry breaking. In fact, several theories beyond the SM, such as Super- 
symmetry and those with extended gauge sectors, require an enlarged Higgs sector in order 
to break the additional symmetries. Here we examine the effect in D^-T^ mixing of models 
with multiple Higgs doublets, with and without flavor conservation, Higgsless models and 
models with scalar leptoquarks. 
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FIG. 11: Left: xd in Vector Leptoquark models as a function of the Vector Leptoquark mass 
MvLQ, with Xpp = 0.1 (soHd line), Xpp = 0.07 (short dash), and Xpp = 0.05 (long dash) for 
P = L or R. The la experimental bounds are as indicated, with the yellow shaded area depicting 
the region that is excluded. 

Right: The present la excluded region in the vector leptoquark mass - coupling parameter plane, 
as well as possible future contours taking xd < (15.0,8.0,5.0,3.0) x 10"^, corresponding to the 
blue dashed, red dashed, cyan dotted, and green dot-dashed curves, respectively. 

A. Flavor Conserving Two-Higgs-Doublet Models 

A simple extension of the SM is to enlarge the Higgs sector by one additional SU{2) 
doublet. We first examine two-Higgs-doublet models that naturally avoid tree- level FCNC 
be requiring that all fermions of a given charge receive their masses from only one Higgs 
doublet |39|]. In one such model, known in the literature as Model II, one doublet (^2) gives 
mass to the up-type quarks, while the down-type quarks and charged leptons receive their 
mass from the other doublet 0i. This is the scenario that is present in supersymmetric 
theories and is in fact required by supersymmetry in order to generate masses for all the 
fermions. Another model, known as Model I, imposes a discrete symmetry such that one 
doublet (^2) generates masses for all fermions and the second (^i) decouples from the fermion 
sector In both cases, each doublet receives a vacuum expectation value Vi, subject to the 
constraint that vf + V2 = v^u- There are five physical scalars in these models, , , A^, 
and H^. The charged Higgs boson can participate in the box diagram for AMp, in exchange 
with the SM Q = —1/3 quarks as shown in Fig. [121 The interactions with the quark 
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FIG. 12: Box diagrams with charged Higgs contributions to D meson mixing. 

sector are governed by the lagrangian 
9 



2^/2Mw 

with Au = cot (3 in both models and 



H [VijUiu^AuUiil - 75)dj + Vijmd^AdUiil + 75)^^] + h.c. 



(75) 



cot /?(tan/?) in Model I(II), where tan/3 = ^2/^1. 



The can have a large contribution 721] to the rare decay h — ^ 57, and in Model II the 
branching fraction for this process sets the bound [73| Mh± ^ 295 GeV at 95% C.L.. This 
constraint is relaxed when other sources of New Physics also contribute to 6 S7, such as 
in supersymmetry. In this case, the lower limit on the charged Higgs mass is 78.6 GeV from 
LEP II data Q. 

It is clear that the contributions to AM^ from the first term in this lagrangian (which 
are proportional to mc,„ cot /5) will only be sizable for extremely small values of tan/5; this 
region is, however, already excluded 74] from, e.g., b ^ S'j and B^-^^ mixing. In addition, 
since the contributions in Model I go as cot (3 multiplied by small mass factors, the effects 
in this case are also restricted to be small. However, the term proportional to mb^s tan f3 in 
Model II has the potential to generate a significant contribution to D mixing in the large 
tan (3 limit. We will thus work in this limit here. Restrictions on the size of tan/3 can be 
obtained by requiring that the ibH^ coupling remain perturbative. If we demand that this 
coupling not exceed the value of the strong coupling constant, gs, we find that tan/3 < 100. 

The effective hamiltonian governing D meson mixing in the large tan/? limit takes the 
form 



GlMl 



2HDM 



W 



47r2 



^ Xi\j{iBxi^ pAnnixi, Xj^xu) + tan^ f3AwH{,Xi, Xj,XH)}Qi 



(76) 



where Aj = VuiV*^ as usual, Xi = m^/M^, the sum extends over i,j = s, b, the loop functions 



can be found in Ref. 



74| and are given in the Appendix. The operator structure is the same 



as in Eq. (!34|) : here, the Qi operator appears due to the presence of the fermion propagator. 
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Note that this structure is quite different than for the case of Bd^s mixing [75| in the large 
tan/3 hmit. This is simply due to the helicity structure of the couplings when the charged 
— 1/3 quarks are internal. The QCD evolution to the charm quark scale is simple and results 
in the factor of ri(mc, Mh±). The resulting contribution to the mass difference is 



X 



D 



X ^ XiXj tan^ IdAuH^Xi, Xj, xh) + tan^ PA^H^Xi, Xj, xh) 

id 



(77) 



Our results for x^^^^^ are displayed in Fig. [13] as a function of tan f3 for various values 
of the charged Higgs mass. We see that the effect is quite small and even at large values 
of tan/3 the contributions from this model are at least an order of magnitude below the 
experimental observation. 
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FIG. 13: x^^^^ in the flavor conserving two-Higgs-doublet model as a function of tan/3 for 
charged Higgs boson masses of m,ii± = 100 , 250, and 500 GeV, corresponding to the solid, dashed 
red, and dashed-dot green curves, respectively. The la experimental bounds are as indicated, with 
the yellow shaded region depicting the region that is excluded. 



In the down quark sector, sizable effects in Bd^s and K meson mixing are obtained in 
the large tan/3 limit from a double penguin contribution with neutral Higgs exchange 3] 
as depicted in Fig. [HI In this limit, an effective Yukawa interaction is induced for the 
down type quarks which includes a contribution from the large Yukawa coupling of the top 
quark. This generates sizable one-loop FCNC couplings for the neutral Higgs in the down 
quark sector. While the same effects occurs in the up quark sector, the term that becomes 
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FIG. 14: The dipenguin diagram with neutral Higgs exchange. 

significant at large tan (3 is now proportional to the down quark Yukawa coupling and hence 
does not generate a sizable effect in D^-T^ mixing 

B. Flavor Changing Neutral Higgs Models 

It is well known that the existence of multiple Higgs doublets can lead in general to tree 
level FCNC transitions [39]. In the down quark sector, there are severe constraints on such 
couplings from Kaon decays, but these do not necessarily lead to equally strong restrictions 
in the up-quark sector. 

The phenomenological requirement that FCNC effects in the down-strange sector must 
be very small can be met in a variety of ways. For example, the imposition of global 



symmetries can make AS = 1 FCNC vanish without affecting the |AC| = 1 sector [76] . 
Another example is the Cheng-Sher ansatz [t^], where the flavor changing couplings of the 
neutral Higgs bosons are given by A/iOj.j^. ~ (v^G^)^/^ym^mjAjj, with the being the 
relevant fermion masses and Aij representing a combination of mixing angles. 

To keep our initial discussion general, we allow for Higgs scalars, which have the 
interactions in the up quark sector 

Cr^X:^^Q^,UR,4>n , (78) 

where Ql represents the left-handed quark doublet and ur is the singlet state. If Mh is the 
mass of the lightest physical Higgs with flavor-changing couplings, the most general effective 
four- fermion hamiltonian just below the Mh scale is 

nn = pGi Qs + G2Q7 + G3 Q,] , (79) 
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where the couphngs Gi,2,3 are model-dependent parameters. As shown in Ref. [78|, 



^1 — X! '^12n^21m^nNA^^ , 
nmN 

nmN 



(80) 



at the Mh scale. Here, AnN refers to the mixing matrix which rotates the Higgs doublets 
$„ to their N neutral physical eigenstates. Matching at the Higgs mass scale relates the 
Wilson coefficients to the three couplings Gi,2,3 via 



H) 



-G 



3 ; 



(81) 



with all other Wilson coefficients being zero. Assuming that Mhj^ > frit for all N and 
computing the evolution of Eq. ( |79l) to fi = rric we obtain 

= 777^ [C3(m,) Q, + C,{m,) Q, + C,{m,) Q, + Cr(m,) Qj + ^(m,) Qg] , (82) 

with 
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[r4(m„ Mh) - r^{m,, Mh)] C^^Mh) , 



r^{m,, Mh) + \ + -^^j '^8("^c, Mh) 



CriM^ 



H) 



[rr{m„ Mh) - rs{m„ Mh)] Cj{Mh) . 



(83) 



8^241 

The Higgs tree-level contribution to xd is found by evaluating the D^-to-D^ matrix 
element, which gives 



X 



(H) _ hflMpBp 



C;{m,) - T] iCim,) + Grim,)) + ^ (G.im,) + Gs{m,)) 
5 5 



4) 



24rBM^ 

where, again, we have used Eqs. (!24l) . (!26l) . Together with Eqs. (IHTj) . (183!) . the above can be 
used to constrain the lightest Higgs mass and associated couplings. As an example, let us 
assume that IGil = \G2\ = at the Higgs mass scale in Eq. (|79ll and Mh is the effective 
mass of the N Higgs scalars. In that case, the restriction on possible values of the effective 
Higgs mass are presented in Fig [151 
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FIG. 15: xd as a function of MujC in models with no natural flavor conservation in the Higgs 
sector. The la experimental bounds are as indicated , with the yellow shaded area depicting the 
region that is excluded. 

We now return to the specific case of the Cheng-Sher ansatz. Here, the neutral Higgs 
bosons can contribute to AM^) through tree-level exchange as well as mediating D meson 
mixing by and t-quark virtual exchange in a box diagram. The restrictions placed on the 
parameter space of this model from the tree-level contribution are computed as described 
above, and are presented for an effective Higgs mass Mh as a function of the coupling 
parameter A^jc in Fig. [161 We see that the form of the couplings, being proportional to the 
light quark masses, result in reduced limits compared to those in Fig. [15] for the general case. 
The box contribution with t-quark exchange is described by the effective hamiltonian 
just below the Mu scale of 

Has = ^^"^"g^.^f + ^a] , (85) 

where x = mf/Mfj, Ftnix) is given in the Appendix, and the vector operators Qis are 
generated from the fermion propagators. The RG evolution and evaluation of the matrix 
elements yields 



The resulting constraints from this contribution are displayed in Fig. [16] in the effective 
Higgs mass - coupling parameter plane. We see that this box contribution only competes 
with those from the tree-level process for large values of Ajj. 
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FIG. 16: Left: la excluded region in the effective neutral Higgs mass - coupling plane for the tree- 
level contribution to D^-D^ mixing in the Cheng-Sher ansatz. Right: la excluded region in the 
effective neutral Higgs mass - coupling plane for the box diagram contribution to D^-D^ mixing in 
the Cheng-Sher ansatz. In both figures, possible future contours taking xd < (15.0 , 8.0 , 5.0 , 3.0) x 
10"^, corresponding to the blue dashed, red dashed, cyan dotted, and green dot-dashed curves, 
respectively, are also displayed. 

C. Scalar Leptoquark Bosons 

Leptoquarks are color triplet particles which couple to a lepton-quark pair and are natu- 
rally present in many theories beyond the SM which relate leptons and quarks at a more fun- 
damental level. Their a priori unknown couplings can be parameterized as X^^/An = Figa. 
Searches for the pair production of scalar leptoquarks at the Tevatron Run II yield the 
bounds j3] 'mLQ> 225 GeV, which are independent of the coupling strength Fi^ 

Scalar leptoquarks participate in D meson mixing via virtual exchange inside a box 
diagram [69;], together with a charged lepton or neutrino. Their interactions are analogous 
to those of R-parity violating supersymmetric models with the terms in the superpotential 
proportional to A'. We thus refer to Section [VIII CI for the details of the analysis for these 
contributions. The resulting constraints on scalar leptoquark are governed by the translation 

= ^'^^^-f^^i^-'^) with = . (87) 
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D. Higgsless Models 



A class of composite Higgs models which has been developed recently to generate a 



naturally light Higgs, employs chiral symmetries of 'theory space' [80| (see also Refs. 40l.l81|). 
Such models involve the 'deconstruction' of higher- dimensional field theories such that the 
low-energy effective field theory resembles the Standard Model but has nice features such as 
the absence of quadratic divergences of the Higgs mass. Here, the Higgs can be interpreted 
as a Goldstone boson of some interaction occuring at higher energies. This approach allows 
for the construction of realistic theories of electroweak symmetry breaking in four spacetime 
dimensions without any higher dimensional interpretation. 

This picture emerges from the AdS/CFT correspondence of the 5-dimensional Higgsless 
model of Csaki et al. ^2\. In the 5-d framework, there is no physical Higgs boson and elec- 
troweak symmetry breaking is generated via the boundary conditions for the 5**^ dimension. 
The gauge symmetry in the higher dimensional space is SU{2)i x SU{2)r x U{1)b~l and 
the right-handed gauge fields receive Planck scale masses. The Kaluza-Klein towers of the 



7 and Z bosons unitarize the WW high energy scattering amplitude [83|] , although there is 
some tension with precision electroweak data as to the precise energy scale that the Kaluza- 
Klein states populate 8j]. The Standard Model fermion fields are localized within the S*'^ 
dimension and also receive their masses from the boundary conditions, with the exact value 
being dependent on their position in the extra dimension 85|. The effects on D mixing 
from this 5-dimensional picture are presented below in Section IVII CI Here, we present our 
results for the AdS/CFT related framework with a composite Higgs. 

The key idea in the composite Higgs picture is that the flavor physics responsible for 



generation of the Yukawa couplings can induce flavor-changing neutral currents 
to charm physics they generically lead to the following effective hamiltonian, 

M2 

^2 



80|. Applied 



H 



W 



C=1,T° 



i^i/ -JJ^^^Llf, C Cl) {ulY C Cl) 



+2 (clsl) 



(44) 



9 
M2 



(ml7^ C Cl) (Mij7^ C cr) 



(88) 



+ (44)' («/?7m C cr) {ur^^ C cr)] 



where g, M are respectively the gauge coupling and gauge boson mass of new flavor gauge 
interactions, and the mixing angles generate different strengths for the gauge coupling. In 



46 



the sum over the color label C, the case C = 1 corresponds to color-singlet interactions, 
whereas C = T'^ refers to color-octet interactions for which T" = A'*/2 are the generators of 
SU (3)f;. In addition, the angles 61 that relate the gauge and mass eigenstates [80;| appear 
in factors oi cl = cos 91 j^, si = sin6'£^, where we take 91 ~ 9c, with 9c being the 
Cabibbo angle. 

The hamiltonian of Eq. (!88l) can be easily transformed to contain the operators from the 
general basis of Eq. f|T6l) , 

= {clslf-^-{cfQ, + CiQ, + C^Q, + CiQ,) , (89) 

where Cf = (3iV, - l)/(2iV,), Cf = rLR{2N, - l)/(2iV,), = -vlr, Cf = ri^Cf , and 
f^LR = (crsI^) / {cIsD . Performing the RG running, we obtain the effective hamiltonian at 
the scale rric with the Wilson coefficients 

Ci(m,) =ri(me,M)Ci(M) , 
C2(m,) =r2(m„M)C2(M) , 

C3(m,) = ^ [r2(m„ M) - rs{m,, M)] C2(M) + r,{m,, M)C^{M) , 

C6(m,) =r6(m„M)C6(M) . (90) 

This, in turn, implies for the mixing amplitude. 



+ ^C^3(m,)(l + 6r/) 



2 (Ci(me) + C6(m,)) + -C2(me) (^- + ^ 



(91) 



where r] is as in Eq. (126|) . The available experimental data can be used to constrain the mass 
of the gauge boson M for different values of the coupling constant as shown in Fig. [TTl 
It is clear that unless g is very small, it is unlikely that any of the gauge bosons of the new 
flavor interactions of Higgsless models will be directly seen at the LHC. 

VII. EXTRA SPACE DIMENSIONS 

Recent speculation that the geometry of spacetime could resolve the hierarchy problem 
have led to theories with extra spatial dimensions that have veriflable consequences at the 



TeV scale. There are several such models 



86| and the size and geometry of the additional 
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FIG. 17: xd as a function of the new gauge boson mass M in Higgless models for g = 0.1 (solid 
line), g = 0.3 (dash-dot line), and g = 0.5 (dashed line). The la experimental bounds are as 
indicated, with the yellow shaded area depicting the region that is excluded. 

spatial dimensions, as well as the field content that is allowed to propagate within them, 
varies between the different scenarios. When the extra dimensions are compactified, the fields 
that reside in the higher dimensional space (known as the bulk) expand into Kaluza-Klein 
(KK) towers of states. The masses of these KK states correspond to the extra dimensional 
components of the bulk field momentum and are related to the bulk geometry. The obser- 
vation, either directly or by indirect effects, of these KK states signals the existence of extra 
dimensions. The properties of the KK states reveal the underlying geometry of the higher 
dimensional spacetime. 

The extra dimensional theories which yield contributions to AF = 2 processes are those 
in which the SM fermion fields reside in the bulk. Here, we consider three such scenarios: (i) 
universal extra dimensions, (ii) localized fermions in a fiat extra dimension, and (iii) warped 
extra dimensions. 




A. Universal Extra Dimensions 



T 

ries 



le possibility of TeV ^-sized extra dimensions naturally arises in braneworld theo- 



87|. By themselves, they do not allow for a reformulation of the hierarchy problem, 
but they may be incorporated into a larger structure in which this problem is solved. The 
scenario which places all Standard Model fields in the bulk is known as Universal Extra 
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Dimensions [88|]. The simplest model of this type contains a single extra dimension com- 
pactified on an S1/Z2 orbifold. Since branes are not present in this case, translational 
invariance in the higher dimensional space would be preserved without the presence of the 
orbifolding. This leads to the tree-level conservation of the extra dimensional momentum of 
the bulk fields, which in turn implies that KK number is conserved at tree-level while KK 
parity, (—1)" where n denotes the KK level, is conserved to all orders in interactions involv- 
ing the KK states. Two immediate consequences of KK number and parity conservation are 
that the KK states must be produced in pairs, and the lightest KK particle is stable and is 



a dark matter candidate 



89| . The former results in a substantial reduction of the sensitivity 



to such states in precision electroweak and collider data. The present bound from Run I at 



90|. 



the Tevatron on the mass of the first KK excitation is of order 250 GeV 

Since all SM fields reside in the bulk in this model, every SM field expands into a KK 
tower of states. The KK reduction of the 5-dimensional fermion fields leaves a chiral zero- 
mode and a vector-like tower of KK states for each flavor. There is one KK tower for each SM 
gauge boson, as well as the Higgs, and additional towers of KK scalars, a^^^ and a^^, which 
correspond to the physical eigenstates of the mixing between the KK towers associated with 
the SM Goldstone fields and the , remnants from the electroweak gauge KK reduction. 
This mixing also generates scalar KK towers which behave as Goldstone fields, which are 
eaten by the gauge boson KK towers and provide masses for the gauge KK states. The 
additional physical scalar KK towers a^^^ do not have zero-modes. The masses of the KK 
states are roughly degenerate and are given at tree-level by 

mn = {ml + n''/Rlfl^ , (92) 

where Rc represents the compactification radius of the extra dimension and mo is the zero- 
mode mass. The KK states clearly become more degenerate with increasing KK-level (in- 
creasing n). These masses are modified jgij] by loop-induced localized kinetic terms and 
non-local radiative corrections. Given that the effect of these mass corrections occurs at 
two- loop order in D meson mixing, one would expect them to have a small effect. We thus 
neglect them in our initial analysis, but will return to this issue at the end of this Section. 

The contributions to D mixing in this model are box diagrams with the boson KK 
tower, its associated KK Goldstone modes G^^, and the a^^, all in exchange with the KK 
towers associated with the d-, s-, and 6-quarks; the zero-mode analogues of these diagrams 
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are shown in Fig. [121 Note that the conservation of KK parity restricts the KK levels of 
the KK quark and boson being exchanged. In addition, only the quark KK towers that are 
even under the Zi symmetry (and thus have a zero-mode) couple to the external zero-mode 
quarks in the box diagram. The relevant hamiltonian at the compactification scale is then 

HuED = E E \.\S{xt\x^;^)Q, , (93) 

n=l i,j 

which has the same structure as that occurring in Eq. ( IM]) . Here, i,j run over (i, s, h. Note 
that at all KK levels, the CKM structure is the same as that in the SM. Using unitarity of 
the CKM matrix, the function S{xf'\x'j^^) becomes 



b[Xi ,Xj ) - 2_^[^PxY{Xi ,Xj )+l^XY{Xa ,Xa j-l'xYyXi ,X^ )-l'xY{Xj ,Xa J 
XY 

(94) 

with xj""* = (m^)^/(m^)^ where now i,j = s,b, and the sum extends over the bosons 



X,Y = W^^, G'^^Y a^^y The functions FxYix^f^i^Y*) ^-re given in the Appendix of Ref. 92|, 



with the appropriate substitutions of quark flavors relevant for D mixing. After the RG 
evolution of the hamiltonian to the charm-quark scale, this leads to 

]\/f2 oo 

^iUED) ^ "fZ^flMnBnnim,, m,) E E ^^\S{xt\ xf) . (95) 

n=l 1,3 

Looking at the expression in Eq. (!92l) for the KK masses, we see that the d- and s-quark 
KK towers are degenerate and the mass splittings between the 6— and rf— , s-quark towers 
are non-zero, yet small, for the flrst couple of KK levels and then effectively vanish for higher 
KK excitations. The GIM cancellation is thus exact in the case of the s-quark KK tower 
contributions, level by level in the KK tower, and leaves a tiny contribution from the flrst few 
6-quark KK states. However, factoring in that ~ (9(10~^), we see that the contributions 
to D mixing from the 6-quark KK states are numerically negligible. Hence, this model is 
not probed by D^-D^ mixing. 

We now return to the case where mass splittings are generated for the KK states via 
localized boundary terms or loop-induced gauge interactions. Since the above one-loop 
contributions to D mixing essentially vanish due to the degeneracy of the KK towers, perhaps 
a non-negligible effect is obtained once the KK degeneracy is lifted. Examining the latter 
effect flrst, we see from Ref. 9l|] that the non-local radiative corrections yield two classes of 



mass splittings for the fermion flelds: {i) a term which is dependent on the gauge couplings 
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FIG. 18: Tree-level gauge KK exchange that mediates neutral meson oscillations. 

and is flavor independent, and (ii) a term which depends on the fermion's Yukawa couphngs. 
The latter term takes the form 

Smf-mf (^m^j , (96) 

where X = (2,4) for fermion (singlets, doublets), respectively, hf is the fermion Yukawa 
coupling, A represents a cut-off scale which absorbs the logarithmic divergences and fi is the 
renormalization scale. If 1/Rc is of order a few hundred GeV, the third generation quark 
doublet and top quark singlet thus receives a correction from the Yukawa term of order 10- 
20 GeV for the first KK state, while the 6-quark singlet KK excitation remains essentially 
unaffected. Given the small CKM factor for the 6-quark KK contributions to D meson 
mixing, and the effectiveness of the GIM mechanism, we find that this mass splitting is not 
enough to generate a sizable contribution to xd. The second possibility of including the 
localized boundary terms holds the promise of inducing large mass splittings between the 
KK states associated with the various quark fiavors. However, these boundary terms may 
take on essentially any value with no predictivity, leaving a virtual continuum of possible 
contributions to D^-D^ mixing. 



B. Split Fermion Models 

In this scenario, the Standard Model fermions are localized at specific points, tji, where 
< < Rc, in extra TeV~^-sized fiat dimensions. The fermions have narrow Gaussian-like 
wavefunctions in the extra dimensions with the width of their wavefunction being much 
smaller than the compactification radius Rc of the additional dimensions. The placement 
of the different fermions at distinct locations in the additional dimensions, along with the 



narrowness of their wavefunctions, can then naturally suppress 93j operators mediating 



dangerous processes such as proton decay and also provide a mechanism for generating the 
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fermion mass hierarchy 94 1. 



This spht fermion scenario is capable of generating large flavor changing neutral cur- 
rents 



95 



96[ |. In contrast to the fermion sector, the gauge bosons are free to propagate 
throughout the extra dimensions. The gauge KK states have cosine profiles which have dif- 
ferent heights at the various distinct fermion locations, generating non-universal couplings 
to different fermion species. This leads to tree-level FCNC as depicted in Fig. [181 with the 
gluon KK states clearly giving the largest contributions. 

With one extra dimension, the coupling of the n*^ KK gluon to a quark localized at the 
scaled position Uq is determined by the overlap of wavefunctions in the additional dimension 



dyi){y)i){y)G^''\y) ^ I dy cos(n7r|/)e-(^-^^)'^'/'^' ~ cos ( e-"''^'/^' , (97) 
Jo \ Rc J 

where a represents the width of the quark's localized wavefunction with a/Rc <^ 1, and yq 
has been normalized to Rc so that it lies in the range < ?/q < 1. The interaction lagrangian 
in the quark mass eigenstate basis is then 

oo 

C = Y: [V2g,G^^-^ {uLYT^V^^d^^V^^u^ + un^T^V^&^W^^un + (u ^ d)] , (98) 

n=l 

where the product V^'^V^ is the usual CKM matrix, the diagonal matrices C{"]j are the 
wavefunction overlaps given above in Eq. ( 1971) . and the factor of \/2 arises from the renor- 
malization of the KK gauge kinetic terms to the canonical value. Ui refers to the set (u,, Q.tj) . 

The effective hamiltonian mediating D meson mixing is given by (taking the contributions 
from the first two generations to be dominant) 

9 CO 1 

'Hs.ut = E 1^ + 20];i2)g2 + 0];[2)Q6) , (99) 

n=l ri 

where t/"^"'' = V^^^C^'^^Vj^ with i = L,R and M„ is the mass of the ra*^ gluon KK state with 
Mn = n/Rc- For the case of one additional dimension, the sum over the gluon KK tower 
converges, and for the scenario with numbers of extra dimensions > 1, the sum is naturally 
cut-off from the finite width of the fermion wavefunction. Performing this sum |96| and 
making use of the unitarity properties of the V^^, we can write the effective hamiltonian at 
the compactification scale as 

ns,iu = \glRl {\Vl,,V}:i,\^F{y^,,yc,)Q, + 2\Vi:,,Vrr2V^^,^^^^^ 

+ \V^iiV^\^?F{y^^,yc^)Q,) , (100) 
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with yu^,ci being the positions of the up- and charm-quark fields, and 



7r2 



Fix,y) = y|x - y\ , 
-7r2 

G{xi,yi,X2,y2) = ^^{\xi-X2\ + \yi-y2\-\xi-y2\-\x2-yi\) ■ (101) 

Although some cancellations could occur between the Qi^ and Q2 terms by finely-tuning 
the quark positions, and thus decreasing the KK gluon contribution to D meson mixing, 
operator mixing from QCD renormalization would spoil this possibility. The RG running of 
the above effective hamiltonian is the same as that performed for the case of flavor changing 
Z' bosons in Section IV Al with the appropriate replacement of the Wilson coefficients, since 
the same operator basis of Qi,2,6 is present. 

In order to explore the magnitude of the KK gluon FCNC effects we examine a single term 
in the above hamiltonian. This will reduce the number of parameters in the computation 
without significantly changing the results. Choosing the term proportional to Qi yields an 
effective hamiltonian at the charm scale of 

-Hspiu = \ n{m,.M)\V}^,,Vr,2?Qi ■ (102) 

Here, Ay = — is the separation between the localized ul and cl quarks, scaled to 
the compactification radius. This leads to a contribution to xd of 

^is,Ht) ^ J_g2Ry^y ^^^^^^ Mm,,V}:i2?fDMnB, . (103) 
yi D 

Figure [19] shows the range of values for clS 3i function of the separation between 

the Ul and cl states for various values of the compactification scale, where = 1/Rc- In 
our numerical work we have used the natural assumption that (Vl)^^ = iVcKM)ij- We see 
that Xp^^**^ vanishes as the separation of the 2 fermions tends to zero as expected. However, 
for most of the range of Ay, we find that compactification scales of order 100 — 500 TeV are 
excluded by the observation of D^-D^ mixing and hence D mixing provides severe constraints 
on the localization of the up-type fermions within this model. Note that these constraints 
are dependent on the choice of values for the elements of the quark diagonalization matrices 
Vl,r, which are a priori unknown, and could be reduced if quark mixing is tiny in the up- 
quark sector. The worst case scenario would be if the are diagonal and all quark mixing 
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FIG. 19: Left: xd in the Split Fermion model as a function of the separation between the left- 
handed u- and c-quark states in the extra dimension for various values of the compactification 
scale. The la experimental bounds are as indicated, with the yellow shaded region depicting the 
region that is excluded. 

Right: la excluded region in the ul — cl separation and compactification scale parameter plane, 
as well as possible future contours taking xd < (15.0,8.0,5.0,3.0) x 10^^, corresponding to the 
blue dashed, red dashed, cyan dotted, and green dot-dashed curves, respectively. 

occurs in the down-quark sector.^ In this case, strong bounds on the compactification scale, 



similar to those presented here, would be obtained from K meson mixing 



95,^ 



C. Warped Geometries 



In the simplest scenario with warped extra dimensions [97|, known as the Randall- 
Sundrum (RS) model, the hierarchy between the electroweak and Planck scales is generated 
geometrically via a large curvature of a single extra dimension. The geometry is that of a 
5-dimensional Anti-de-Sitter space (AdSs), where the extent of the 5*^ dimension is y = vrrc 

is the compactification radius), and every slice of the additional dimension corresponds 
to a 4-d Minkowski metric. Two 3-branes reside at the boundaries of the AdSs slice, with 
the 3-brane located at the fixed point y = ttTc being known as the TeV brane, while the 



^ This is frequently the case in models of quark mass matrices where the up-quark mass matrix is taken to 
be diagonal and all mixing is assigned to the down-quark sector. A rationale for this is given, for example. 



in Ref. [38 1. 
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opposite brane at the other boundary y = is referred to as the Planck brane. Within this 
framework, gravity is locahzed about the Planck brane, and electroweak symmetry breaking 
can take place either with the Higgs field being localized on or near the TeV brane, or via 
boundary conditions imposed at the fixed points as in the Higgsless models discussed above. 
The FCNC effects considered here are independent of this choice. 
FCNC effects are induced 198, 



99| when the SM fermions and gauge bosons are local- 



ized in the warped 5 dimension 



100 



101 



102| . As in the case of flat TeV ^-sized extra 



dimensions with split fermions discussed above, the observed fermion masses and mixings 
are automatically explained by the geometry, with the 5-dimensional Yukawa couplings all 
being of order unity. Localizing the light fermions near the Planck brane results in small 
4-dimensional Yukawa couplings for these fields, whereas if the top-quark field is localized 
near the TeV brane a large 4-d top Yukawa coupling is induced. This localization scheme 
also naturally suppresses higher dimensional flavor changing operators that are problematic 
when the SM is confined to the TeV brane. This fiavor breaking fermion localization leads 
to FCNC interactions via non-universal couplings of the zero-mode fermions to the gauge 
boson KK states. Since the SM gauge bosons are localized near the TeV brane (in order 
to acquire their masses) and have exponentially decaying wavefunctions towards the Planck 
brane, we expect FCNC in the light quark sector to be suppressed. 
The action for fermion fields in the RS bulk is given by 101 1 



S = j (fx dyVO Q^7^Dm^ + sgn{y)Mf^^ + h.c^ , (104) 

where G represents the determinant of the 5-dimensional metric. Dm is the covariant deriva- 
tive in curved space, and = V^'^'j'^ with V^^ being the inverse vierbein. The parameter 
of importance to us here is Mf which is the 5-dimensional bulk mass for the fermion /. It 
given by Mf = kcf, where k is the parameter describing the curvature of the AdSs space and 
is of order of the 5-d Planck scale. The constants c/ indicate the position of the fermion's 
localized wavefunction in the bulk, with cj > 1/2 (c/<l/2) corresponding to the fermion 
being localized near the Planck (TeV) brane. These constants determine the fiavor structure 
of the theory. 

The KK decomposition of the bulk fermion fields yields the normalized zero-mode wave- 
function (a discrete symmetry ensures that the zero-mode fields are chiral), 



/^°^ = ^ ^cYfC-^f^y . (105) 
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Range of cj 


1 

2 

1 

27rkrc 


1 

C/ < 2 - « 

> i + e 



TABLE IV: The asymptotic behavior of the square of the parameter Yf for various localization 
points {cf} of the fermion's wavefunction. 

The asymptotic behavior of the Yf on the locahzation parameters Cf are hsted in Table HVl 
We see that these factors become exponentially small when the fermions are localized near 
the Planck brane. In the basis where the 5-d bulk masses, Mf, are diagonal the fermion 
Higgs interactions yield the 4-dimensional Yukawa couplings, 

for the zero-mode fermions in terms of the 5-d Yukawa couplings A5. We take the elements 
of A5 to be complex and of order unity. Note that kr^ ~ 11.3 in order to resolve the gauge 
hierarchy problem. The elements of the matrices that diagonalize the up and down quark 
fields to their 4-d mass eigenstates have magnitude 



Y, 



fhi 



Y, 



Wr 



CKM\ij 



(107) 



with L —>■ R for the matrices that diagonalize the right-handed fields. For the elements with 
j < i one should interchange i ^ j. 

The wavefunctions for the gauge KK states are given by the first order Bessel functions 
Ji , Yi, and the mass of the n*^ g auge KK mode is M„ = Xnke'^^^" where x„ is related to 



the roots of Bessel functions 



lOOj]. The first few values of a;„ are 2.45, 5.57, 8.70, and 1 1.84 



Prec ision electroweak data places severe bounds on the masses of the gauge KK states |100l . 
103!]. However these bounds can be reduced to Mi > 3 TeV if the gauge symmetry i n th e 
bulk is expanded to SU{2)i x SU(2)r x U(1)b-l^ which restores custodial symmetry |l04 |. 
The couplings of these states to the zero- mode fermions, {Cj-"^}, are determined by the 
overlap of their wavefunctions in the additional dimension. They are given (as a ratio to the 
SM coupling) by 



9SM 



2TikrcYflf^ , 



(108) 
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where /I"'' is an integral over Ji Bessel functions, and are given explicitly in Ref. 102l | and in 



the Appendix. As displayed in Fig. [201 this coupling weakens substantially as the gauge KK 
level, n, increases for Cf < 1/2, while for > 1/2 the couplings tend to a small fixed value 
for all KK levels. The interaction lagrangian for the gluon KK states in the quark mass 
eigenstate basis is as given in Eq. ( l98l) with the substitution of the prefactor \/2 ^Jl-akrc^ 
which arises from the renormalization of the KK gauge kinetic terms to the canonical value. 



BO 




FIG. 20: The coupling strength, scaled to the SM strong coupling constant, of the zero-mode 
fermions to the first five gluon KK excitations (as labeled) as a function of the fermion bulk mass 
parameter c/. 



D^-D^ mixing is then mediated via tree-level flavor changing interactions of the KK 
gauge boson states as depicted in Fig. [TSl In analogy to the previous section, the effective 
hamiltonian for this process is given by (for the exchange of gluon KK states which yield 
the largest contribution) 

2Tlkrc 2 \ ^ 1 (jju(n)Uju{n) ^ r^j Mn)\ j ju{n) ^ j ju{n)\ j ju{n) ^ \ (TnQ\ 
riRS - -^9s 2^ ITp ['^L{cu)'^Liuc)Ql + ^'^Licu)'^R{uc)Q2 + 1^ Ricu)'~^ R{uc)Q6 ) ' (1^9) 



where Ul%^ = V^'^^^CfVl^. Writing this explicitly for Ul"^'"' yields 



n=l n 



u{n) 



u 



u(n) 



L{uc) - ^L{uj)'-^jk ^L{kc)"jk - ^L(uj)^ fj^L{jc) 



(110) 



smce the df^ are diagonal, and where the index j sums over the generations. Looking at 
the asymptotic values of Yf in Table HVt we see that f/^l"]) ^^^y sizable when the fermion 
is localized towards the TeV brane. Unitarity of the V^j^ results in 



n 



RS 



"g^, {C,{M^)Q, + C2(M„)g2 + CeiM^)Qe] 



;iii) 
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where Mi is the mass of the first gluon KK excitation, and with the Wilson coefficients 
being given by 

n=l n 

+ <i2)^r(22)K(/i?)^ - Yiiii-^r]) Wi3)V^i^(32)K(/ff )^ - Y^iii-jr] 
+ <i2)V«"(22)K^.(4?)^ - Yiiii-jr]) , 

CO i\r2 

C,{M,) = 2.kr, E ^ «3)^i^(32)K?.(/£V - 

n=l n 

+ <12)V^«(22)K',(/S^)' - , (112) 

with Ml/M"^ = xl/x'^ where the Bessel function roots described above. 

The RG evolution proceeds as in Section III Al and results in the effective hamiltonian at 
the charm quark scale 

Urs = ^ (Ci(m,)Qi + C2K)Q2 + C3K)Q3 + C^{ni,)Q^) , (113) 

where additional operators have been generated due to mixing in the RG evolution. The 
evolved Wilson coefficients at the charm scale are as given in Eq. (l52ll with the appropriate 
substitution of Mz' Mi. Upon evaluating the matrix elements we obtain the contribution 
to Xd from warped extra dimensions 

= Jl^^^^ (^i^iK) + CeK)] - ^C2K) + ^C3K)) (114) 

in the modified vacuum saturation approximation. Here, we have taken the factor rj of 
Eq. ([26]) to be unity. 

To obtain numerical results, we need to specify the fermion locations in the warped 
dimension. We examine three popular scenarios in the literature that correctly generate the 
4-d Yukawa hierarchy for the SM fermions. As mentioned above, localizing the fields near 
the UV (ultraviolet or Planck) brane generates an exponentially small 4-d Yukawa coupling. 
In all three models, all of the light quarks are localized such that their bulk mass parameters 
take on values with Cf > 1/2. Special attention must be paid to the localization of the third 
generation quarks; in order to generate a large top-quark mass, the corresponding SU(2) 
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Model I 


Model II 


Model III 




> 1/2 


0.6 


0.5 




> 1/2 


0.6 


1.4 




> 1/2 


0.52 


0.5 


CCR 


> 1/2 


0.52 


0.53 


CtL 


0.45 


0.4 


0.46 


CtR 





0.3 


on the IR brane 



TABLE V: The values of the bulk mass parameters for the three models described in the text. 



singlet field is usually taken to reside close to the TeV brane. The third generation SU(2) 
doublet fields and 6-quark SU(2) singlet field are located as close to the IR (infrared, or TeV) 
brane while maintaining consistency with the experimentally determined Zbb coupling. The 
three scenarios that we follow are fairly uniform in their treatment of the third generation, 
differing only slightly in the location of the SU(2) top-quar k sin glet. The scenarios are: 
(I) A study of flavor physics in the Randall- Sundrum model jl05 |. (II) A scenario that has 
been constructed in order to generate fermion masses within the 5-d picture of Higg sless 



models 85j, (III) The up-quark singlet field is taken to lie even closer to the UV brane 106 1 
in order to solve the strong CP problem with warped geometries. The numerical values of 
the bulk mass parameters are summarized in Table |V] for the three cases. 

Our results for x^^ for these three models are presented in Fig. EH where as above, we 
assume the quark diagonalization matrices take on CKM-like values. In the figure, the dot- 
dashed green, dashed red, and solid black curves correspond to the bulk mass parameters of 
Model I, II, and III, respectively. We see from Fig. |20] that fermions localized towards the 
Planck brane have very small couplings to the KK gluon states and thus do not substantially 
contribute to x^^ . This simplifies the expressions in Eq. f lll2p in this only the 

tL,R terms have sizable contributions. Looking at the figure we see the mass of the first 
gluon KK excitation is constrained to lie ^ 1 — 2 TeV, which is essentially the same value 
as the bound obtain ed f rom the precision electroweak data in warped models with bulk 
custodial symmetry [l04| . Lastly, we recall from discussion in the previous section, that 
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these constraints can be evaded if the matrix which diagonahzes the up-quark sector is 
essentially diagonal. 




FIG. 21: The contribution to xd from a warped extra dimension with the SM fields in the bulk as a 
function of the mass for the first gluon KK excitation. The la experimental bounds are as indicated, 
with the yellow shaded region depicting the region that is excluded. The curves correspond to Model 
I (dashed-dot, green), Model II (dashed, red) and Model III (solid) as described in the text. 



VIII. EXTRA SYMMETRIES 



In this section, we focus on supersymmetry. Weak scale supersymmetry is a possible 
solution to the gauge hierarchy problem, leads to the unification of the gauge couplings at 
high energies, and provides a natural Dark Matter candidate. It is thus a very well motivated 
theory of physics beyond the SM. Supersymmetry is an extension of the Poincare symmetry, 
relating fermions and bosons at a fundamental level. All SM particles have supersymmetric 
partners sparticles') with the same mass and gauge interactions, but with spin differing by 
one-half unit. Since the supersymmetric particles have yet to be discovered, we know that 
supersymmetry is broken; in this section, we will be agnostic as to which supersymmetry 
breaking mechanism Nature may have chosen. We note that in non-broken supersymmetry 
the rates for all loop-induced processes would vanish due to an exact cancellation between 
the SM and supersymmetric contributions. It is thus due to the breaking of supersymmetry 
that contributions to FCNC are generated in these theories. 

Here, we examine the contributions to D^-D^ mixing in four supersymmetric scenar- 
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ios: the Minimal Supersymmetric Standard Model (MSSM), models with alignment in the 
quark-squark mass matrices, models with R-parity violating couplings, and Split Supersym- 
metry. Other scenarios with extended non-supersymmetric symmetries have been considered 
elsewhere in this paper. 



A. Minimal Supersymmetric Standard Model 

As the name implies, the MSSM is the simplest version of supersymmetry as it contains 
the minimal number of new particles. The SM fermions are placed in chiral supermultiplets, 
the SM gauge bosons lie in vector supermultiplets, and the Higgs sector takes the form of 
the flavor conserving two-Higgs-doublet Model II discussed above. A discrete symmetry, 
R-parity, is imposed to forbid unwanted terms in the superpotential that would mediate 
proton decay at a dangerous level. Conservation of R-parity implies that only pairs of 
sparticles can be produced or exchanged in loops. Collider searches for direct squark and 
gluino pair production place the bound rriq^g > 330 GeV [l^ in the MSSM with gravity 
mediated Supersymmetry breaking. 

As mentioned above, we will not assume any particular supersymmetry breaking mecha- 
nism in our discussion, and so we employ a model independent parameterization of all possi- 
ble soft supersymmetry breaking terms. This soft supersymmetry breaking sector generally 
includes three gaugino masses, trilinear scalar interactions, as well as Higgs and sfermion 
masses, and thus contains many potential sources of flavor violation. Here, we are interested 
in the flavor violating sources that arise in the up-squark sector. In what is known as the 
super-CKM basis, the squark fields are rotated by the same matrices that diagonalize the 
quark masses, giving rise to non-diagonal squark mass matrices. The squark propagators are 
then expanded suc h that th e non - diag onal mass terms result in mass insertions that change 



the squark flavor |l07l . llOSl . Il09l . Ill0||. This source of flavor violation differs from that of 
the SM and many NP models discussed earlier. Here, the quark-squark-gaugino neutral 
couplings are flavor conserving, while flavor violation arises from the non-diagonality of the 
squark mass propagators. The 6x6 mass matrix for the Q = +2/3 squarks can be divided 
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FIG. 22: Contributions to D mixing from mass insertions in the squark propagator in MSSM. 
N ,M ,P, and Q label the hehcity {L,R). 



into 3x3 sub-matrices, 



ML ML 



\MUM^ 



(115) 



RR j 



and the mass insertions can be parameterized in a model independent fashion as 



m-, 



(116) 



Here, z,j are flavor indices, M,N refers to the helicity choices LL, Li?, RR, and irtq rep- 
resents the average squark mass. Although this source of flavor violation is present in 
general, and in particular in models with gravity mediated supersymmetry breaking, it can 
be avoided if supersymmetry is broken by gauge or anomaly mediation. These mass inser- 
tions are thought to be small in the MSSM, but can be large in non-minimal supersymmetric 
models. 

In this scenario, the virtual exchange of squarks and gluinos in the box diagrams depicted 
in Fig. [22] can have a strong contribution to D^-D^ mixing. Note that the second diagram 
in the figure is due to the Majorana nature of the gluino. The effective hamiltonian at the 
supersymmetric scale is given by 

„2 8 



n 



a: 



MSSM 



2m? 

'^'"'q 1=1 



;ii7) 



where all eight operators in the ind epen dent basis contribute. The matching conditions at 



the supersymmetric mass scale are 



ml 



1 ,..,,2 



Ci(m|) = ^ (5^2)^ [4^/i(^) + ll/2(^)] , 



1^ 
1 



C2{m'q) = m2)LR iSi2)RL 15/2(a:) - (5^,),, (6^,)^^, [2x/i(x) + 10/2(x)]} 

_ 1 
"1^ - 9 



CsimI) = ^ m,),, (6-,)^^ [42xA(x) - 6Ux)] - (6-,)^^ (^^2)^^ HUx)} 
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C4K-) = -(5r2)L37x/i(x) 

C5K-) = l(5r2)k ^/l(^), 



(118) 



K-) = Y^W2)L [4X/1(X) + 11/2(X)] 



C7K-) = Y^W2)L37x/i(x), 
where x = m?/m?, with rrig being the mass of the gluino. The functions fi{x) and 



are given in the Appendix. Note that these conditions are symmetric under the interchange 
L ^ R. We also note that the NLO expressions for these matching conditions have been 



computed in Ref. 112[ |. 

The RG evolution to the charm-quark scale results in 



Ci(mc) = ri(mc,mg)Ci(mg) 



C^{mc) 
Crime) 



1) ' 



1) ' 



Cz(m^ = - [r2(mc, rriq) - r^{mc, m^)] C2{mq) + rsirric, mq)C3{mq) , 



V241 



[r5(mc, niq) - r^{mc, nig)] 



15 

Ciirrig) + —C^irrig) 



+ 2 rUg) + r^irric, rUg)] C^{mg) 



1 



[r4(mc, ruq) - r^iruc, rUg)] [Ci{mg) + GAC^img)] 



8V24T 

+ - [r-iimc, rrig) + r^im^, nig)] C^{mg) , 



rQ{mc,mg)Ce{mg) , 
g 

-j== [rsim^, nig) - r-jiym^, nig)] 



15 



Ciiniq) + —Cs{mg) 



+ - [rjinic, nig) + r^{nic, nig)] Cjinig) , 



[r7(mc, mg) - r^{nic, m,-)] [C7(mg) + 64C8(mg)] 



8V24T 

+ - [r7(mc, m^) + r^{nic, nig)] Cs{nig) , 



which agrees in form with that in Ref. 



(119) 



113l |. Here, we have assumed that the squarks and 



gluinos are integrated out at roughly the same scale. Upon evaluating the matrix elements 
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in the modified vacuum saturation approximation we obtain the MSSM contribution to xd, 



^{MSSM) _ 0!s fo^D^D 



2 '"2 5 7 



q ^ D 



(120) 



6 

Here, we have taken the factor rj of Eq. ( 126|) to be unity. 

Our results for presented in Figs. [231 1211 [23 In these figures, we show 

contours for the absolute value of the up-charm squark mass insertions for various helicities as 
a function of the ratio mg/rriq for different average squark masses. These contours correspond 
to xd = (11.7,15.0,3.0) X lO^'^ in the three figures. In figure [23| the region above the 
contours represents the current la excluded region. In these figures we take one or two of 
the mass insertions to be non-vanishing, as indicated. Due to the L <-> symmetry of the 
matching conditions, the constraints on |5"2Ui |5i2|r-R5 well as |5i2Ui? l^ul RL 
are identical. We see that D meson mixing restricts the up and charm squark masses to be 
degenerate at the 1 — 10% level for most of the parameter space. We n ote that our results 



1141 ]. 



numerically agree with those recently computed by Ciuchini et al. 

There are several other contributions to D meson mixing within the MSSM. These are 
all mediated via box diagrams with internal sparticle exchange and we now discuss each 
one in turn, (i) The exchange of any of the 4 neutralinos (^ = 1)4) with the up and 
charm squarks. This contribution proceeds via mass insertions in the squark propagators 
with flavor diagonal quark-squark-neutralino couplings as in the case of internal squark- 
gluino exchange discussed above. Since the couplings are of weak interaction strength in 
this case, the magnitude of this contribution is suppressed by the ratio g^/g^ compared to 
the squark-gluino results and is thus numerically insignificant, (ii) The exchange of one of 
the neutralinos and one gluino with the up and charm squarks. This again proceeds via the 
non-diagonal squark mass insertions at a rate of g"^ / g1 compared to the pure gluino-squark 
contribution. Although larger than the pure neutralino-squark contribution, it is still a 
sub-leading effect, (iii) The exchange of charginos xf = 1; 2) and all three down-type 
squarks inside the box diagram. Here, the squark propagators are diagonal (mass insertions 
do not contribute in this case since the internal squarks are Q = —1/3) and the flavor 
violation is given by the CKM structure of the quark-squark-chargino vertices. However, 
the Q = —1/3 squarks are constrained to be highly degenerate from their contributions (with 
gluino exchange) to K , , Bg meson mixing. Thus a supersymmetric-GIM mechanism is in 
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FIG. 23: The constraints on the absolute value of the mass insertions with different helicities as a 
function of the mass ratio rrig/mq for various values of the average squark mass. The la excluded 
region, corresponding to xd < 11.7 x 10^'^, lies above the curves. 

effect, yielding nearly exact cancellations, and rendering this contribution negligible. This 
is in contrast to the chargino-squark contributions to K ,B(i ,Bs meson mixing, where the 
potentially non-degenerate stop squark participates and can induce large contributions, (iv) 
The charged Higgs contribution of the two-Higgs-Doublet model of type II discussed in 
Section IVI A[ As shown in that section, these contributions are numerically small, even in 
the case of large tan /3. In summary, we see that all other supersymmetric contributions to 
D^-D^ mixing are numerically insignificant compared to the squark-gluino exchange. It is 
interesting to note that stop-squarks do not contribute to D meson mixing. 
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FIG. 24: Contours, corresponding to xd = 15.0 x 10^'^, for the absolute value of the mass insertions 
with different helicities as a function of the mass ratio rrig/niq for various values of the average 
squark mass. The region above the curves corresponds to larger values of xd. 

B. Quark-Squark Alignment Models 



As we saw in the previous Section, in the MSSM there is a new "flavor problem," namely, 
how to keep the contributions from the supersymmetric particles to FCNC as small as the 
observations. The conventional solution is to impose constraints, such as those derived 
above, of (i) degeneracy in the squark sector (except for the special case of stop squarks), 
i.e. the diagonal sub-matrices Mil and Mrji in Eq. f lll6p should be proportional to the 
unit matrix, and (ii) the non-diagonal sub-matrices M^/j should be proportional to the 
corresponding quar k ma trix. 

Nir and Seiberg 115l | have proposed an alternative to this picture where the quark and 
squark mass matrices are approximately aligned with each other. Their proposal is as follows: 
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FIG. 25: Contours, corresponding to xd = 3.0 x 10~^, for the absolute value of the mass insertions 
with different helicities as a function of the mass ratio mg/mq for various values of the average 
squark mass. The region above the curves corresponds to larger values of xd. 

if for some symmetry reason the matrices corresponding to the squark mass insertions, 5mni 
are themselves diagonal, then the squark contributions to FCNC vanish, regardless of the 
mass spectrum of the squarks. Corrections to this approximation are expected to remain 
tolerably small and it should be possible to simultaneously diagonalize the quark mass 
matrices and the squark mass-squared matrices while essentially preserving flavor diagonal 
gluino interactions. 

Within this framework, it is somewhat problematic to satisfy the constraints from 
mixing. Specific implementations of this proposal, based on Abelian horizontal symmetries, 
restrict the supersymmetric contributions to Kao n mi xing via a unique structure for the 
down quark mass matrix using holomorphic zeros 116l |. This implies that Cabibbo mixing 
between the first and second generation quarks must be induced by mixing in the up-quark 



67 



sector, which in turn leads to sizable supersymmetric contributions to D^-D^ mixing. In 
this case, mixing in the up-charm squark sector gives 

{Sll)^.= ^ ^- ^^e^ , (121) 

where 6c is the Cabibbo angle, while the {6lr)uc mass insertions can naturally remain small. 
Mirroring the above discussion for MSSM, this leads to the effective Hamiltonian that me- 
diates D mixing 



= 77^CiK-)Qi , (122) 



with 



2ml 



C^{m,) = ^i6-,,)U4xf^{x) + nUx)] , (123) 



where fi^2ix) with x = m^/m?- are again given in the Appendix. The RG evolution is simple 
and yields 

x\^' = —^ ri(m„mg)(:7i(mg) . (124) 

The bounds on i^^\2}LL from the current measurement of D meson mixing are given in the 
upper left-hand panel of Fig. [231 Using Eq. (11211) above, this results in the constraint on 
squark and gluino masses of ( assu ming rrig for simplicity) g > 2 TeV, which agrees 



with the results in Refs. 



12 



114l |. This would exclude early discovery of Supersymmetry 



at the LHC, but leaves a discovery window with higher luminosities as the LHC detectors 
are expected to have a search reach of g up to 2.5 — 3.0 TeV with 300 fb^^ of integrated 
luminosity. 



C. Supersymmetry with R-Parity Violation 

The conventional gauge symmetries of supersymmetry allow for the existence of additional 
terms in the superpotential that violate baryon and lepton number. The assumption of R- 
parity conservation in the MSSM prohibits these terms, ensuring that baryon and lepton 
number are conserved, and forbids related dangerous operators, e.g., those that mediate 



proton decay. However, it is possible to construct alternative discrete symmetries |117l |. 
such as baryon-parity or lepton-parity, that allow terms which violate either baryon or 
lepton number, but not both. These symmetries also forbid unwanted operators, and there 
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is no strong theoretical motivation to prefer R-parity over these alternative scenarios. The 
R-parity violating terms in the superpotential can be written as 

Wr^ = ^XijkLiLjEk + X[,kLiQjDk + ^Kj,UiD,Dk . (125) 

i,j, k are generation indices and symmetry demands i j {j 7^ k) in the terms proportional 
to A (A"). The quantities L, E, Q, D, U in Eq. ( I125P are the chiral superfields in the MSSM, 
and the SU (2) l, SU (3)c indices have been suppressed. A bilinear term may also be present, 
but it can be rotated away and will not be considered here. The lepton number violating 
terms, A and A', cannot exist simultaneously with the AB 7^ term containing A". The 
A' terms have the same structure as the couplings for scalar leptoquarks, as discussed in 
Section IVI CI This model still contains the minimal superfield content, but leads to a 
markedly different supersymmetric phenomenology as sparticles can now be produced singly 
and can mediate FCNC at tree-level. 

The superfields are in the weak basis and should be rotated to their mass eigenstates. 



The AL 7^ A' term becomes 118l | 



WR^ = ~K^,[N,V,iDi~E,Uj]Dk , (126) 

with the definition 

Kjk = Krs^rj^lk ■ (127) 

Here, and are the matrices which rotate the left-handed up- and right-handed down- 
quark fields to their mass basis. Written in terms of component fields, the second term in 
this superpotential contains the interactions 

W^v = a:,, {V,i [vldld}, + dl4z/l + (4)*(^1)^4] (128) 
-el4<-«i4el-(4)*(el)X} , 

where the second line involving the up-quark sector is relevant for D^-D^ mixing. 

Constraints on the size of these R-parity violating couplings have been ob taine d in the 



literature. These limits are derived from considerations of various processes 1191] such as 
charged current universality, semi-leptonic meson decays, rare meson decays, atomic parity 
violation, double nucleon decay, neutron oscillations, and Z boson decays. A compilation of 
the 2a bounds on the couplings relevant for D^-D^ mixing are given in Table IVTI In addition. 
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TABLE VI: 2(T constraints on the R-parity violating couplings which participate in D mixing. 
Here, A; = 1,2,3 with the exception that k ^ j, where j represents the middle index, for the A" 
couplings. All numbers are scalec by the factor (mj^ ^/lOO GeV). Details of the derivation of these 
restrictions are given in Refs. 3) 113] • 



the recently improved upper bou nd on the branching fraction for the process D 



1201 1 yields the stringent restriction [7'] on the product of 

2 



of i3 < 7.4 ■ 10"^ from CLEO-c 
couplings ~X[,,X[,, < 0.003 (m,-^^/(100 GeV) 

It is possible for the quark flavor rotations to generate flavor violation in the down- or 
up-quark sectors, but not both. In the case where the flavor rotations occur in the up-quark 
sector only, large flavor changing effects are expected in the D meson s ystem and the limits 



on the R-parity violating couplings shown in Table IVll become modified [1191]. However, this 
scenario is rather model dependent, we will adopt a more conservative, model-independent 
formalism in the following. 

For the lepton number violating coupling A', the first and third terms in the second line 
of Eq. (128) mediate D^-D^ mixing via box diagrams where either the pair [ii^i — dn^k) or 
{^T,,i — dj i.k) are exchanged internally with the assignment of the generational index j = 



1, 2 1181 ]. The corresponding Feyman diagrams are depicted in Fig. [261 Note that there are 
no tree-level contributions as in the case of meson mixing in the down-quark sector. This is 
described at the high mass scale by the effective hamiltonian 



n 



Rr, 



1287r' 



'XK2kK 



ilk) 



1 1 



7X1% 



m-, 



Q 



1 ) 



(129) 

where the dependence on the operator Qi is induced due to the fermion propagator. This 
interaction will yield constraints on the product of couplings A^a/c-^iifc- Here, we have assumed 
that only one set of the R-parity violating couplings Xi2kKik {i-^-, only one value of i and 
k) is large and dominant. This is equivalent to saying that, e.g., both sleptons and both 
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FIG. 26: Contributions to D^-D^ mixing from the A' superpotential terms in supersymmetric 
models with R-parity violation. 

down-type quarks being exchanged in the first box diagram shown in Fig. [22] are from the 
same generation. In general, this need not be the case and, for example, the coupling factor 
would then be the product ^'i2k'^'mik^'m2nKin^ with, e.g., the set of iL,i,dR^k,(^L,m,dR^n being 
exchanged. 

Matching at the SUSY scale yields the Wilson coefficient 

^iK-) = ^(AkA:j^(l + J^l ■ (130) 



Computing the evolution to the charm-quark scale yields 



= ^^C^{m,)Q, , (131) 

with 

Ci{mc) = ri{mc,mg)Ci{mq) . (132) 

Evaluating the appropriate matrix element gives the D mixing contribution for the R-parity 
violating A' terms, 

d.R 

Taking m^^ . ~ mj^ ^ for simplicity, we obtain the constraint 

iK2kKlkf < (cxpt) QGtt^Fj, 

which yields numerically 

A:,X, < 0.085v^ (^^^) . (135) 

We find that relaxing our assumption on the slepton mass and taking m^^ m^^^ strength- 
ens this bound at most by a factor of 3.7 when m^-^ . = 100 GeV. In computing the RGE 
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evolution we used the value rriq = 500 GeV and find little sensitivity in the evolution on the 



squark mass once it is above the current experimental limit from HERA of ~ 300 GeV |13 |. 
It is trivial to scale our result in Eq. (11351) to compare to the limits in Table I VI I which are 
based on setting mi =100 GeV. Taking, mr = 500 GeV, we see that the bounds on 
K2kKik from D^-D^ mixing are a factor of 50 (250) times stronger than those in the Table 
for i = 2 (i = 3). 

The full numerical results for xd and the constraints obtained in the R-parity violating 
coupling, squark mass parameter plane are presented in Fig. |27] in the limit m^^ ~ m ^ . 
We see that D meson mixing provides stringent constraints on R-parity violating couplings. 
These bounds can be directly translated to constraints on the couplings of scalar leptoquarks 
as discussed in a previous Section. 
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FIG. 27: Left: xd in supersymmetry with R-parity violation as a function of the product of R- 
parity violating couplings )^i2k^'iik taking "t-j^^ = "i^^ with "ij^ = 300,500,1000 and 2000 
GeV corresponding to the solid, green dashed, red dotted, and blue dashed-dot curves, respectively. 
The 1(7 experimental bounds are as indicated, with the yellow shaded region depicting the region 
that is excluded. Right: lex excluded region in the R-parity violating coupling - squark mass plane, 
as well as possible future contours taking xd < (15.0,8.0,5.0,3.0) x 10~^, corresponding to the 
blue dashed, red dashed, cyan dotted, and green dot-dashed curves, respectively. 

The baryon number violating A" couplings can also contribute to D^-D^ mixing. We 
remind the reader that they cannot exist simultaneously with the lepton number violating 
terms in the superpotential. They participate in D mixing via dR-qnaik and dn exchange 
in the box diagram. The formalism is analogous to the AL ^ case above. The effective 
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hamiltonian at the SUSY mass scale is 

1 



2 ' 2 

ml m- 



Recall that symmetry dictates j ^ k. After completing the RGE evolution as described 



Qi . (136) 



above, and assuming the charged —1/3 squarks are degenerate, we have 



X 



° 967r2m| To 



ri(me,mj^)(A"^.fcA2^.J . (137) 



This yields the constraint 

KA<^-^^^^'{-^^) , (138) 

which mirrors that for the lepton number violating scenario. 
D. Split Super symmetry 



Lastly, for completeness, we briefly discuss the case of Split Supersymmetry [l2l[ |. This 
scenario postulates that supersymmetry breaking occurs at a very high scale, irts ^ 1000 
TeV. The scalar particles all acquire masses at this high scale, except for a single neutral 
Higgs boson, whose mass is either finely-tuned or is preserved by some other mechanism. 
Split Supersymmetry proponents argue that this tuning may, indeed, be present in Nature, 
perhaps being related to the cosmological constant problem (which suffers an even greater 
degree of fine-tuning). The fermions in this theory, including the gauginos, are assumed 
to be protected by chiral symmetries and thus can have weak-scale masses. This feature 
preserves the gauge coupling unification found in supersymmetric models, and provides a 
natural Dark Matter candidate in the lightest neutralino. One important consequence of this 
scenario is that since all the scalar fields are present at only a very high scale, they decouple 
from physics at the TeV scale and their contributions to FCNC in the flavor sector are 
negligible. Since all contributions to D^-D^ mixing in supersymmetry involve the internal 
exchange of scalar quarks, we expect these effects to essentially vanish in this scenario. 

IX. CONCLUSIONS 

The recent BaBar and Belle findings on D^-D^ mixing have brought the long standing 
search for this phenomenon to a successful conclusion, although much remains to be done. 
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Compared to mixing in the other flavor sectors, the observed value for charm (c/Eq. (jl])) is 
by far the smallest, 

xk ^ 0.47 , Xb, ^ 0.776 , xb, ^ 26. , xb ^ 0.009 . (139) 

In our opinion, the measured value for x-d is in accord with expectations of the Standard 
Model, with the proviso that hadronic (rather than quark-level) effects are the dominating 
influence (c/Sect. III). We have argued that the relatively small magnitude of charm mixing 
could afford New Physics an enhanced chance to compete successfully with the Standard 
Model. It should be kept in mind, however, that the current experimental value of xd is 
relatively imprecise and that the SM theoretical determination contains hadronic uncer- 
tainties. These facts tend to frustrate the attempt to disentangle any potentially large NP 
contribution from that of the SM. 

By design, our study of NP contributions has addressed a rather broad spectrum of 
possibilities. We have avoided playing favorites among the NP models contained in this 
paper, letting the results speak for themselves. Since the average reader is unlikely to be 
conversant with the details of such a large array of NP models, our presentation has been 
pedagogical in nature. We have tried to precede any formula for xj^^'' with a summary of 
the relevant background. 

A work such as this is meant to constrain the parameter spaces of NP models. The case 
of mixing is especially interesting because the intermediate states which generate the 
D^-io-D^ transitions are distinct from those occurring in K, Bd and mixing. Typically, 
NP parameters will involve the masses of yet-to-be-discovered particles and their coupling 
strengths to ordinary matter. In some cases (Left-Right Symmetric Model, Split Supersym- 
metry. Universal Extra Dimensions, Flavor Conserving Two-Higgs Doublets) we have found 
that the NP model will not generate a D^-D^ signal at the observed level for any values of its 
parameters. More often, however, this is not the case and for some models (Split Fermions, 
Flavor Changing Neutral Higgs) the constraints can be strong. 

The main quantitative conclusions for this work appear in the set of flgures which appear 
throughout the paper. For convenience, we have compiled a summary of our results in 
Table IVnl using the la value xd < 11.7 ■ 10~^ to mark the boundary between allowed and 
excluded regions. Such a list is by nature approximate, and we refer the reader to the body 
of the paper for a more precise presentation of our results. 
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iVlodei 


Approximate Constraint 


Fourth Generation (Fig. ^ 


iKb'Kfe'l • W < 0.5 (GeV) 


g = -1/3 Singlet Quark (Fig. HI) 


S2-ms < 0.27 (GeV) 


Q = +2/3 Singlet Quark (Fig. El) 


|A„e| < 2.4-10-4 


Little Higgs 


Tree: See entry for Q = —1/3 Singlet Quark 




Box: Region of parameter space can reach observed xd 


Generic Z' (Fig. ED 


Mz'/C > 2.2 • 10^ TeV 


Family Symmetries (Fig. [8|) 


mi// > 1.2 • 10^ TeV (with mi/m2 = 0.5) 


Left-Right Symmetric (Fig. [9]) 


No constraint 


Alternate Left-Right Symmetric (Fig. [TO]) 


Mr > 1.2 TeV {ruD^ = 0.5 TeV) 




{Am/mD,)/MR > 0.4 TeV-^ 


Vector Leptoquark Bosons (Fig. [TT]) 


MvLQ > 55(App/0.1) TeV 


Flavor Conserving Two-Higgs-Doublet (Fig. [T3l) 


No constraint 


Flavor Changing Neutral Higgs (Fig. [T5]) 


mn/C > 2.4-103 TeV 


FC Neutral Higgs (Cheng-Sher ansatz) (Fig. [TB]) 


mH/\Auc\ > 600 GeV 


Scalar Leptoquark Bosons 


See entry for RPV SUSY 


Higgsless (Fig. [TT]) 


M > 100 TeV 


Universal Extra Dimensions 


No constraint 


Split Fermion (Fig. \T9\i 


M/\Ay\ > (6 - 10^ GeV) 


Warped Geometries (Fig. [2T|l 


Ml > 3.5 TeV 


Minimal Supersymmetric Standard (Fig. [23l) 


|(5]^2)lr,rl| < 3.5 - 10-2 for m ~ 1 TeV 




|(53'2)ll,rr| < -25 for m ~ 1 TeV 


Supersymmetric Alignment 


m > 2 TeV 


Supersymmetry with RPV (Fig. [27]) 


A'i2feA'iifc/m,-^^ < 1.8 - 10-3/100 GeV 


Split Supersymmetry 


No constraint 
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TABLE VH: Approximate constraints on NP models from mixing. 



We recommend further experimental study of this subject on two fronts. First, of course, 
is the need to reduce error bars in the measured values of yr, and especially xb- Equally 
important is continuing the search for evidence of CP violation in mixing for the system. 
CP violation provides an interesting contrast with mixing because it provides an inde- 
pendent arena for competition between the SM and NP signals. There is especially room 
for improvement in the SM analysis of charm CP violation, and work on this is underway. 
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APPENDIX A: COLLECTED FORMULAE 



Here, we collect the formulae used throughout the manuscript to compute the contribu- 
tions to D^-D^ mixing in the various New Physics models. 
Inami-Lim: (from Ref. 32j) 

The loop functions, first calculated by Inami Lim, apply for several NP scenarios discussed 
in the text. For a contribution from two internal quarks of the same flavor in the box diagram, 
the loop function is 



S{x) 



X 



9 



1 

4 ^ 4(1 -x) 



2 1 



x] 



Sx'^ 



2 1 



\nx 



X] 



(Al) 
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and for two quarks of different flavors, 



Inx,- 



i3j 2 J 



1 3 _ 3 

4^2(1-Xi) ~4(1-Xi)2 

^ 1 



(A2) 



Lztt/e Higgs: (from Ref. [47|) 

Tfie loop function for the case where the mirror fermions and heavy gauge bosons are 
exchanged in the box diagram is given by 



1 



4 



{Zi - Zj){l - ZiY 



z] log Zj 



{Z, - Zj){l - Zj^ 

3 / 1 



l-2zi + 



Zj^Zj 



+ 



zf log Zi 



Zj log Zj 



Zi){l- Zj) {z, - Zj){l - Zi^ {zi - Zj){l - ZjY 
3 / 1 



(A3) 



7T + 



Z■Z^\0gZ■ 



z'jZj log Z'j 



100a \{1 - zl){l - z'j) {z - I - Zj){l - zl)^) {zi - Zj){l - z'j)\ 

_3_ / logQ zf log Zj 

10 [{a - 1)(1 - - z^) {z, - Zj){l - z,)(l - 

Zj log Zj \ 



{Z, - Zj){l - Zj){l - z'j) ^ 



with 



M,, 



m 



M,, 



tan^ 9,, 



(A4) 



Left-Right Symmetric Model: 

The loop function with one Wl and one Wr boson being exchanged in the box diagram 



IS 



J(x,/?) 



x(3\n(3 



X + X In X 



4l-/5)(l-/?x)2 (l-x)(l-/3x) ■ 



(A5) 



Charged Higgs: (from Ref. |74| ) 

The loop functions with one and one SM W boson and with two being exchanged 
in the box diagram are 



A}in{x,y) 



X 



2xy , X 
^ In 



AwR{x,y) = 2x^ 



X + y 

{x - {x - y)^ y_ 

1 y\ny 



(x^ — Inx 

{x-y){l-x) ' (x-?/)2(l-y) ' (x-|/)2(l-x) 



+ 



+ 



X 



+ 



y^ \ny 



+ 



x(x + xy — 2y) Inx' 



4 \{x — y){l — x) {1 — y){x — y)'^ (x — 2/)^(l — x) 



• (A6) 
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Cheng-Sher box: 

The loop function in the Cheng-Sher ansatz with flavor changing neutral Higgs bosons 
for a top quark and neutral Higgs being exchanged in a box diagram is 

-1 



FtH(x) 



1 — X (1 — xy 

Universal Extra Dimensions: 



Inx — 4x + 3 + 2 In X 

+ - 



2(1 -x¥ 



(AT) 



The expressions for the case of Univers al Ex tra Dimenions are given in full in Ref . 92 1 . 
Warped Extra Dimensions: (from Ref. 102l |) 

The coupling of two zero-mode fermions to the n*^ gauge boson KK state, A^'^\ relative 
to the SM coupling strength is 



(jffA 



QSM 



2Txkrr 



I -2c 



f 



]^ ^-TTkrc{l-2cf ) 



' ci.z(--/)#^4^^^#^ , (A8) 



where the roots x„ for the gauge boson KK spectrum are given by 



Jl{Xn) + XnJ[iXn) + [YliXn) + XnYKx^)] = 



and an is defined by 



Ji{mn/k) + {mn/k)J[{mn/k) 
run/k) + {mn/k)Y{{mn/k) 



(A9) 



(AlO) 



Supersymmetry: (from Ref. 

The loop functions from squark and gluino exchange in a box diagram with squark mass 
insertions are given by 

6(1 + 3x) \nx + x^ - - 9x + 17 



/i(^) 
/2(x) 



6(1 

6a;(l + x)lnx — x^ — 9x^ + 9x + 1 
3(l-a;)5 



(All) 
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